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ABSTRACT 

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined 
with the WMAP large-angle polarization, constrain the scalar spectral index to = 0.9603 ± 0.0073, ruling out exact scale invariance at over 5cr. 
Planck establishes an upper bound on the tensor-to-scalar ratio at r < 0.1 1 (95% CL). The Planck data shrink the space of allowed standard infla- 
tionary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential 
models of degree « > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, 
obtaining d«s/dlnA: = -0.0134 ± 0.0090. We verify these conclusions through a numerical analysis, which makes no slow-roll approximation, 
and carry out a Bayesian parameter estimation and model-selection analysis for a number of inflationary models including monomial, natural, 
and hilltop potentials. For each model, we present the Planck constraints on the parameters of the potential and explore several possibilities for 
the post-inflationary entropy generation epoch, thus obtaining nontrivial data-driven constraints. We also present a direct reconstruction of the 
observable range of the inflaton potential. Unless a quartic term is allowed in the potential, we find results consistent with second-order slow-roll 
predictions. We investigate whether the primordial power spectrum contains any features. A penalized likelihood approach suggests a feature 
near the highest wavenumbers probed by Planck at an estimated significance of ~ 3cr after correction for the "look elsewhere" effect. We find 
that models with a parameterized oscillatory feature improve the fit by ^xls ~ 10; however, Bayesian evidence does not prefer these models. 
We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with Planck bounds on /ml. 
Planck constrains with unprecedented accuracy the amplitude and possible correlation of non-decaying isocurvature fluctuations. The fractional 
primordial contributions of CDM isocurvature modes of the types expected in the curvaton and axion scenarios have upper bounds of 0.25% and 
3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelated isocurvature component 
can improve the xlg by ~ 4 due to a moderate tension between / < 40 and higher multipoles. Nonetheless, the data are consistent with adiabatic 
initial conditions. 
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1. Introduction 

This paper, one of a set associated with the 2013 release of data 
from the Planck^ mission (Planck Collaboration 1, 2013), de- 
scribes the implications of the Planck measurement of cosmic 
microwave background (CMB) anisotropics for cosmic infla- 
tion. In this first release only the Planck temperature data re- 
sulting from the 'nominal' mission is used, which includes 2.6 
surveys of the sky. The interpretation of the CMB polarization 
as seen by Planck will be in a later series of publications. This 
paper exploits the data presented in Planck Collaboration XV 
(2013); Planck Collaboration XVll (2013). Other closely related 
papers are the estimates of cosmological parameters in Planck 
Collaboration XVI (2013) and investigations of non-Gaussianity 
in Planck Collaboration XXIV (2013). 

In the early 1980s inflationary cosmology, which postu- 
lates an epoch of nearly exponential expansion, was proposed 
in order to resolve a number of puzzles of standard Big Bang 
cosmology such as the entropy, flatness, and monopole prob- 
lems (Starobinsky, 1980; Guth, 1981; Albrecht & Steinhardt, 
1982; Linde, 1982, 1983). During inflation, cosmological fluc- 
tuations resulting from quantum fluctuations are generated and 
can be calculated using the semiclassical theory of quantum 
fields in curved spacetime (Mukhanov & Chibisov, 1981, 1982; 
Hawking, 1982; Guth & Pi, 1982; Starobinsky, 1982; Bardeen 
etal., 1983; Mukhanov, 1985). 

Cosmological observations pie-Planck are consistent with 
the simplest models of inflation within the slow-roll paradigm. 
Recent observations of the CMB anisotropics (Story et al., 2012; 
Bennett et al., 2012; Hinshaw et al., 2012a; Hou et al., 2012; 
Das et al., 2013) and of large-scale structure (Beutler et al., 
201 1; Padmanabhan et al., 2012; Anderson et al., 2012) indicate 
that our Universe is very close to spatially flat and has primor- 
dial density fluctuations that are nearly Gaussian and adiabatic 
and are described by a nearly scale-invariant power spectrum. 
Pre-Planck CMB observations also indicate that the amplitude 
of gravitational waves, with a nearly scale-invariant spectrum 
generated in the early Universe (Grishchuk, 1975; Starobinsky, 
1979; Rubakov et al., 1982; Fabbri & Pollock, 1983), is small. 

Most of the results in this paper are based on the two- 
point statistics of the CMB as measured by Planck, exploit- 
ing the data presented in Planck Collaboration XV (2013); 
Planck Collaboration XVI (2013) and Planck Collaboration 
XVll (2013). The Planck results testing the Gaussianity of the 
primordial CMB component are described in the companion pa- 
pers Planck Collaboration XXlll (2013); Planck Collaboration 
XXIV (2013); Planck Collaboration XXV (2013). Planck finds 
values for the non-Gaussian /nl parameter connected to the 
CMB bispectrum that are consistent with the Gaussian hypoth- 
esis (Planck Collaboration XXIV, 2013), and this result has im- 
portant implications for inflation. The simplest slow-roll infla- 
tionary models predict a level of bispectral non-Gaussianity /nl 
of the order of the slow-roll parameters and therefore too small 
to be detected by Planck. 

The paper is organized as follows. Section 2 reviews in- 
flationary theory emphasizing in particular those aspects used 
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later in the paper. In Section 3 the statistical methodology and 
the Planck and other astrophysical likelihoods are described. 
Section 4 studies constraints on slow-roll inflation and studies 
their robustness under generalizations of the minimal assump- 
tions of our baseline cosmological model. In Section 5 Bayesian 
model comparison of several inflationary models taking into ac- 
count the uncertainty from the end of inflation to the beginning 
of the radiation dominated era is carried out. Section 6 recon- 
structs the inflationary potential over the range corresponding 
to the scales observable in the CMB. In Section 7 a penalized 
likelihood reconstruction of the primordial perturbation spec- 
trum is performed. Section 8 reports on a parametric search for 
oscillations and features in the primordial scalar power spec- 
trum. Section 9 examines constraints on non-canonical single- 
field models of inflation including the /nl measurements from 
Planck Collaboration XXIV (2013). In Section 10 constraints on 
isocurvature modes are established, thus testing the hypothesis 
that initial conditions were solely adiabatic. We summarize our 
conclusions in Section 1 1 . Definitions of the most relevant sym- 
bols used in this paper can be found in Tables 1 and 2. 

2. Lightning review of inflation 

Before describing cosmic inflation, which was developed in the 
early 1980s, it is useful to review the state of theory prior to 
its introduction. Lifshitz & Khalatnikov (1963) first wrote down 
the equations for the evolution of linearized perturbations about 
a homogeneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) spacetime within the framework of general rel- 
ativity. The general framework adopted was based on two as- 
sumptions: 

(i) the cosmological perturbations can be described by a single 
component fluid, at least at very early times; 

(ii) the initial cosmological perturbations were statistically ho- 
mogeneous and isotropic, and also Gaussian. 

These are the simplest — but by no means unique — assumptions 
for defining a stochastic process for the initial conditions. 
Assumption (i), where only a single adiabatic mode is excited, is 
just the simplest possibility; in Sect. 10 we shall describe isocur- 
vature perturbations, where other available modes are excited, 
and report on the constraints established by Planck. Assumption 
(ii) is a priori more dubious given the understanding at the time. 
An appeal can be made to the fact that any physics at weak cou- 
pling could explain (ii), but at the time these assumptions were 
somewhat ad hoc. 

Even with the strong assumptions (i) and (ii), comparisons 
with observations cannot be made without further restrictions on 
the functional form of the primordial power spectrum of large- 
scale curvature inhomogeneities K, Vnik) oc k"'^'^, where n^ 
is the (scalar) spectral index. The notion of a scale-invariant 
(i.e., Hs - 1) primordial power spectrum was introduced by 
Harrison (1970); Zeldovich (1972); and Peebles & Yu (1970) 
to address this problem. These authors showed that the scale- 
invariant power law was consistent with the crude constraints on 
large- and small-scale perturbations available at the time; how- 
ever, other than its mathematical simplicity, no compelling the- 
oretical explanation for this ansatz was put forth. An important 
current question, addressed in Sect. 4, is whether n^-l (i.e., ex- 
act scale invariance) is consistent with the data, or whether there 
is convincing evidence for small deviations from exact scale in- 
variance. Although the inflationary potential can be tuned to ob- 
tain Hs = 1, inflationary models generically predict deviations 
from «s = 1, usually on the red side (i.e, «s < 1). 
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Table 1. Cosmological parameter definitions 



Parameter Definition 



fit Baryon fraction today (compared to critical density) 

Q.C Cold dark matter fraction today (compared to critical density) 

Omc Approximation to the angular size of sound horizon at last scattering 

T Thomson scattering optical depth of reionized intergalactic medium 

A'eff Number of massive and massless neutrinos 

Sum of neutrino masses 

Yp Fraction of baryonic mass in primordial helium 

Q.f: Spatial curvature parameter 

Wde Dark energy equation of state parameter (i.e., p/p) (assumed constant) 

"R Curvature perturbation 

I Isocurvature perturbation 

Px = k^l^kl"^ /'Ztt^ ■ Power spectrum of X 

Ax X power spectrum amplitude (at k, - 0.05 Mpc"') 

Ms Scalar spectrum spectral index (at k, = 0.05 Mpc"') 

dWs/dlnA; Running of scalar spectral index (at k, = 0.05 Mpc"') 

(fiis/d Ink^ Running of running of scalar spectral index (at k, = 0.05 Mpc"' ) 

r Tensor-to-scalar power ratio (at k, = 0.05 Mpc"') 

n, Tensor spectrum spectral index (at k, = 0.05 Mpc"') 

dnjdlnk Runrung of tensor spectral index (at k, = 0.05 Mpc ') 



2. 1. Cosmic inflation 

Inflation was developed in a series of papers by Starobinsky 
(1980); Guth (1981); Albrecht & Steinhardt (1982); and Linde 
(1982, 1983). (See also Sato (1981) and Brout et al. (1978) for 
interesting precursors.) By generating an equation of state with 
a negative pressure (i.e., w = p/p ~ -1) before the radia- 
tion epoch, inflation solves a number of cosmological conun- 
drums (the monopole, horizon, smoothness, and entropy prob- 
lems), which had plagued all cosmological models extrapolat- 
ing a matter-radiation equation of state all the way back to the 
singularity. Such an equation of state with a large negative pres- 
sure and the resulting nearly-exponential expansion are obtained 
from a scalar field, the inflaton, with a canonical kinetic term, 
slowly rolling in the framework of Einstein gravity. 

The homogeneous evolution of the inflaton field, ^, is gov- 
erned by the equation of motion 



defined as 



'4>{t) + 3//(f)0(f) + = 0, 

and the Friedmann equation 



1 /I 



+ y{4>) 



(1) 



(2) 



Here H - a /a is the Hubble parameter, the subscript ^ denotes 
the derivative with respect to 0, Mpi = (8;7rG)"'^^ is the reduced 
Planck mass, and V is the potential. The evolution during the 
stage of quasi-exponential expansion, when the scalar field rolls 
slowly down the potential, can be approximated by neglecting 
the second time deiivative in Eq. 1 and the kinetic energy term 
in Eq. 2, so that 



V(<f>) 



3Mrf2 



(3) 
(4) 



Necessary conditions for the slow-roll described above are ey <K 
1 and |?7v| <k 1, where the slow -roll parameters ey and rjy are 



M\VI 
pi 

2y2 ' 
y 



(5) 
(6) 



An analogous hierarchy of parameters, called the Hubble 
flow function (HFF) slow-roll parameters, measures instead the 
deviation from an exact exponential expansion. This hierarchy 
is defined as e\ - -H/H^, e,+i = eiliHei), with i > 1. By using 
Eq. 3, we have that » ey, €2 ~ -2t]v + A-ey 

2.2. Quantum generation of fluctuations 

Without quantum fluctuations, inflationary theory would fail. 
Classically, any initial spatial curvature or gradients in the scalar 
field, as well any inhomogeneities in other fields, would rapidly 
decay away during the quasi-exponential expansion. The result- 
ing Universe would be too homogeneous and isotropic compared 
with observations. Quantum fluctuations must exist in order to 
satisfy the uncertainty relations that follow from the canonical 
commutation relations of quantum field theory. The quantum 
fluctuations in the inflaton, and in the transverse and traceless 
paits of the metric, are amplified by the nearly exponential ex- 
pansion, yielding the scalar and tensor primordial power spectra, 
respectively. 

Many essentially equivalent approaches to quantizing the 
linearized cosmological fluctuations can be found in the origi- 
nal literature (see, e.g., Mukhanov & Chibisov, 1981; Hawking, 
1982; Guth & Pi, 1982; Starobinsky, 1982; Bardeen et al., 1983). 
A simple formalism, which we shall follow here, was intro- 
duced by Mukhanov (1988); Mukhanov et al. (1992) and Sasaki 
(1986). In this approach a gauge-invariant inflaton fluctuation Q 
is constructed and canonically quantized. This gauge-invariant 
variable Q is the inflaton fluctuation 6(p(t, x) in the uniform cur- 
vature gauge. The mode function of the inflaton fluctuations 
d(p{t, x) obeys the evolution equation 



(a6,Pk)" + {k^ -j\(a6,f>k) = 0. 



(7) 
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Table 2. Conventions and definitions for inflation physics 



Parameter Definition 



^ Inflaton 

Vi(j)) Inflaton potential 

a Scale factor 

/ Cosmic (proper) time 

5X Fluctuation of X 

X = dX/dt Derivative with respect to proper time 

X' = dX/dr] Derivative with respect to conformal time 

X^ = dX/d(f) Partial derivative with respect to if> 

Mpi Reduced Planck mass (= 2.435 x 10'** GeV) 

Q Scalar perturbation variable 

h*'^ Gravitational wave amplitude of (+, x)-polarization component 

X, X evaluated at Hubble exit during inflation of mode with wavenumber k, 

X^ X evaluated at end of inflation 

£v = Mp|V^/2V^ First slow-roll parameter for V((f>) 

Tjv = M^^V^^/V^ Second slow-roll parameter for V((p) 

^2 = M^^V^V^^^/V^ . . Third slow-roll parameter for V{(p) 

rUy = M^^VjV^ip^^ / . Fourth slow-roll parameter for V((f>) 

£i=-H/H- First Hubble hierarchy parameter 

e„+i = e„IHe„ (n + l)st Hubble hierarchy parameter (where n > 1) 

N(t) = f 'dt H Number of e-folds to end of inflation 

Sa Curvature field perturbation 

6s Isocurvature field perturbation 



where z - a^/H. The gauge-invariant field fluctuation is directly 
related to the comoving curvature perturbation- 

K = -H-!- . (8) 

Analogously, gravitational waves are described by the two 
polarization states {+, x) of the transverse and traceless parts of 
the metric fluctuations and are amplified by the expansion of the 
Universe as well (Grishchuk, 1975). The evolution equation for 
their mode function is 

(fl/2;'T + (fc'-^j(«/'r) = o. (9) 

Early discussions of the generation of gravitational waves dur- 
ing inflation include Starobinsky (1979); Rubakov et al. (1982); 
Fabbri & Pollock (1983); Abbott & Wise (1984). 

Because the primordial perturbations are small, of order 
10"^, the linearized equations 7 and 9 provide an accurate de- 
scription for the generation and subsequent evolution of the cos- 
mological perturbations during inflation. In this paper we use 
two approaches for solving for the cosmological perturbations. 
Firstly, we use an approximate treatment based on the slow- 
roll approximation described below. Secondly, we use an almost 
exact approach based on numerical integration of the ordinary 
differential equations 7 and 9 for each value of the co-moving 
wavenumber k. For fixed k the evolution may be divided into 
three epochs: (i) sub-Hubble evolution, (ii) Hubble crossing evo- 
lution, and (iii) super-Hubble evolution. During (i) the wave- 
length is much smaller than the Hubble length, and the mode os- 
cillates as it would in a non-expanding Universe (i.e., Minkowski 
space). Therefore we can proceed with quantization as we would 

2 Another important quantity is the curvature perturbation on uniform 
density hypersurfaces f (in the Newtonian gauge, ^ = —ij/ - Hdp/p, 
where i// is the generalized gravitational potential) which is related to 
the perturbed spatial curvature '•^^R = -4V-^/a^. On large scales ( ~ R. 



in Minkowski space. We quantize by singling out the positive 
frequency solution, as in the Bunch-Davies vacuum (Bunch & 
Davies, 1978). This epoch is the oscillating regime in the WKB 
approximation. In (iii), by contrast, there are two solutions, a 
growing and a decaying mode, and the evolution becomes inde- 
pendent of k. We care only about the growing mode. On scales 
much larger than the Hubble radius (i.e., k <k aH), both cur- 
vature and tensor fluctuations admit solutions constant in time.-' 
All the interesting, or nontrivial, evolution takes place between 
(i) and (iii) — that is, during (ii), a few e-folds before and after 
Hubble crossing, and this is the interval where the numerical in- 
tegration is most useful since the asymptotic expansions are not 
valid in this transition region. Two numerical codes are used in 
this paper, ModeCode (Adams et al., 2001; Peiris et al., 2003; 
Mortonson et al., 2009; Easther & Peiris, 2012), and the inflation 
module of Lesgourgues & VaUcenburg (2007) as implemented in 
CLASS (Lesgourgues, 2011; Bias et al., 2011).** 

It is convenient to expand the power spectra of curvature and 
tensor perturbations as 

/ , + i dn,/dlnjHn(/i;/A:,)+i d-«Jdln/t^(ln(jt//t.))^+- 

W) = aJ-J " (10) 
v] ' ^^^^ 

where As (At) is the scalar (tensor) amplitude and («t), 
drts / din kidnj dink) and d^Wj/dlnfe^ are the scalar (tensor) 
spectral index, the running of the scalar (tensor) spectral index, 
and the running of the running of the scalar spectral index, re- 
spectively. 

^ On large scales, the curvature fluctuation is constant in time when 
non-adiabatic pressure terms are negligible. This condition is typically 
violated in multi-field inflationary models. 

^ http : / /zu server 2 .star.ucl.ac.uk/ ~ hi r any a /ModeCode / , 
http : //class-code . net 



Planck Collaboration: Constraints on inflation 



5 



The parameters of the scalar and tensor power spectra may 
be calculated approximately in the framework of the slow-roll 
approximation by evaluating the following equations at the value 
of the inflation field 0» where the mode kt = a*//* crosses the 
Hubble radius for the first time. (For a nice review of the slow- 
roll approximation, see for example Liddle & Lyth (1993)). The 
number of e-folds before the end of inflation, Nt, at which the 
pivot scale kt exits from the Hubble radius, is 



Nt 



V 



(12) 



where the equality holds in the slow-roU approximation, and 
subscript 'e' refers to the end of inflation. 

The coefficients of Eqs. 10 and 1 1 at their respective leading 
orders in the slow-roll parameters are given by 



As 
A, 

"s - 1 



V 



2V 

pi 

Irjv - 6ev 
-2ev 



d«Jdlnfe K -16evT]v + 24ey + 2^1 
dnjdlnk -Aevrjv + Se^ 
d\jdlnk^ ^ -1924 + I92el7]v - 32evi]l 
- 24£vfv + '^mfv + 2c7v, 



(13) 

(14) 

(15) 
(16) 
(17) 
(18) 

(19) 



where the slow-roll parameters ey and rjy are defined in Eqs. 5 
and 6, and the higher order parameters are defined as follows 



and 



3 _ K^l^' 



y3 



(20) 



(21) 



In single field inflation with a standard kinetic term, as dis- 
cussed here, the tensor spectrum shape is not independent from 
the other parameters. The slow-roll paradigm implies a tensor- 
to-scalar ratio, at the pivot scale, of 



16e * -Sn, 



(22) 



referred to as the consistency relation. This consistency relation 
is also useful to understand how r is connected to the evolution 
of the inflaton: 



1 

— dA^V^. 
V8 Jo 



(23) 



The above relation, called the Lyth bound (Lyth, 1997), im- 
plies that an inflaton variation of the order of the Planck mass 
is needed to produce r > 0.01. Such a threshold is useful to 
classify large and small field inflationary models with respect to 
the Lyth bound. 



2.3. Ending inflation and the epoch of entropy generation 

The greatest uncertainty in calculating the perturbation spectrum 
predicted from a particular inflationary potential arises in estab- 
lishing the correspondence between the comoving wavenumber 
today, and the inflaton energy density when the mode of that 
wavenumber crossed the Hubble radius during inflation (Kinney 
& Riotto, 2006). This correspondence depends both on the infla- 
tionary model and on the cosmological evolution from the end 
of inflation to the present. 

After the slow-roll stage, ip becomes as important as the cos- 
mological damping term 3i/0. Inflation ends gradually as the 
inflaton picks up kinetic energy so that w is no longer slightly 
above -1, but rather far from that value. We may arbitrarily 
deem that inflation ends when w - -1/3 (the value dividing 
the cases of an expanding and a contracting comoving Hubble 
radius), or, equivalently, at ey a; 1, after which the epoch of 
entropy generation starts. Because of couplings to other fields, 
the energy initially in the form of scalar field vacuum energy 
is transferred to the other fields by perturbative decay (reheat- 
ing), possibly preceded by a non-perturbative stage (preheating). 
There is considerable uncertainty about the mechanisms of en- 
tropy generation, or thermalization, which subsequently lead to 
a standard w = 1/3 equation of state for radiation. 

On the other hand, if we want to identify some fc, today with 
the value of the inflaton field at the time this scale left the hori- 
zon, Eq. 12 needs to be matched to an expression that quantifies 
how much kt has shrunk relative to the size of the horizon be- 
tween the end of inflation and the time that mode re-enters the 
horizon. This quantity depends both on the inflationary potential 
and the details of the entropy generation process, and is given by 



A^. X 71.21 - In 
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(24) 



where pend is the energy density at the end of inflation, p± is 
an energy scale by which the Universe has thermalized, oo^o is 
the present horizon scale, Vhor is the potential energy when the 
present horizon scale left the horizon during inflation, and Wint 
characterizes the effective equation of state between the end of 
inflation and the energy scale specified by pth. In predicting the 
primordial power spectra at observable scales for a specific in- 
flaton potential, this uncertainty in the reheating history of the 
Universe becomes relevant and can be taken into account by al- 
lowing Nt to vary over a range of values. Note that vv'int is not 
intended to provide a detailed model for entropy generation, but 
rather to parameterize the uncertainty regarding the expansion 
rate of the Universe during this intermediate era. Nevertheless, 
constraints on Wi„i provide observational limits on the uncertain 
physics during this period. 

The first two terms of Eq. 24 are model independent, with 
the second term being roughly 5 for kt = 0.05 Mpc If ther- 
malization occurs rapidly, or if the reheating stage is close to 
radiation-like, the magnitude of the last term in Eq. 24 is < 1 . 
For most reasonable inflation models, the fourth term is (9(1) and 

the third term 10, motivating the commonly assumed range 

50 < A^t < 60. Nonetheless, more extreme values on both ends 
are in principle possible (Liddle & Leach, 2003). In the figures 
of Sect. 4 we will mark the range 50 < A^* < 60 to guide the 
reader's eye. 
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2.4. Perturbations from cosmic inflation at higfier order 

To calculate the quantum fluctuations generated during cosmic 
inflation, a linearized quantum field theory in a time-dependent 
background can be used. The leading order is the two-point cor- 
relation function 

<iR(ki) iR(k2)) = {2nf ^'P'Rik) 6\U + ^2), (25) 

but the inflaton self-interactions and the nonlinearity of Einstein 
gravity give small higher-order corrections, of which the next- 
to-leading order is the three point function 

<??(ki) 9?(k2) 9?(k3)) = {2nfBK(h,k2, ^3)5^(ki+k2+k3) , (26) 

which is in general non-zero. 

For single-field inflation with a standard kinetic term in 
a smooth potential (with initial fluctuations in the Bunch- 
Davies vacuum), the non-Gaussian contribution to the curva- 
ture perturbation during inflation is 0{ev , Iv) (Acquaviva et al., 
2003; Maldacena, 2003), i.e., at an undetectable level smaller 
than other general relativistic contributions, such as the cross- 
correlation between the integrated Sachs-Wolfe effect and weak 
gravitational lensing of the CMB. For a general scalar field 
Lagrangian, the non-Gaussian contribution can be large enough 
to be accessible to Planck with /nl ~ c^^ (Chen et al., 
2007), where Cs is the sound speed of inflaton fluctuations (see 
Sect. 9). Other higher order kinetic and spatial derivative terms 
contribute to larger non-Gaussianities. For a review of non- 
Gaussianity generated during inflation, see, for example, Bartolo 
et al. (2004a); Chen (2010), and the companion paper Planck 
Collaboration XXIV (2013). 

2.5. Multi-field models of cosmic inflation 

Inflation as described thus far assumes a single scalar field that 
drives and terminates the quasi-exponential expansion, and also 
generates the large-scale curvature perturbations. When there 
is more than one field with an effective mass smaller than H, 
isocurvature perturbations are also generated during inflation, 
by the same mechanism of amplification due to the stretching 
of the spacetime geometry (Axenides et al., 1983; Linde, 1985). 
Models with an component inflaton 0, can be best analysed in 
terms of an average field cr moving along the classical trajectory 
(&^ = Yjf=i 4>'i)^ whose fluctuations are related to the curvature 
perturbation according to 

^ = -j/ ^^='f^' , (27) 

where 2, is the gauge-invariant field fluctuation associated with 
(pj. The - 1 normal directions are connected to - 1 isocur- 
vature perturbations 6sij according to 

6sij = . (28) 

cr 

If the trajectory of the average field is curved in field space, then 
during inflation both curvature and isocurvature fluctuations are 
generated with non-vanishing correlations. 

Isocurvature perturbations can be converted into curvature 
perturbations on large scales, but the opposite does not hold 
(Mollerach, 1990). If such isocurvature perturbations are not 
totally converted into curvature perturbations and persist after 
nucleosynthesis, these can lead to observable effects in CMB 



anisotropics and in structure formation. In Sect. 10, we present 
the Planck constraints on a combination of curvature and isocur- 
vature initial conditions, and the implications for important two- 
field scenarios, such as the curvaton (Lyth & Wands, 2002) and 
axion (Lyth, 1990) models. 

Isocurvature perturbations may lead to a higher level of non- 
Gaussianity compared to a single inflaton with a standard kinetic 
term (Groot Nibbelink & van Tent, 2000). There is no reason to 
expect the inflaton to be a single-component field: the scalar sec- 
tor of the Standard Model and its extensions contain numerous 
scalar fields. 



3. Methodology 

3. 1. Cosmological model and parameters 

The parameters of the models to be estimated in this paper fall 
into three categories: (i) parameters describing the initial pertur- 
bations, i.e., characterizing the particular inflationary scenario in 
question; (ii) parameters determining cosmological evolution at 
late times (z < 10"*); and (iii) parameters that quantify our un- 
certainty about the instrument and foreground contributions to 
the angular power spectrum. These will be described in Section 
3.2.1. 

Unless specified otherwise, we assume that the late-time cos- 
mology is the standard flat six parameter ACDM model, whose 
energy content consists of photons, baryons, cold dark matter, 
neutrinos (assuming A^etr - 3.046 effective species, one of which 
is taken to be massive with a mass of niy - 0.06 eV), and a 
cosmological constant. The primordial helium fraction, Fp, is 
set as a function of Q.^lv' and A^eif according to the Big Bang 
nucleosynthesis consistency condition (Ichikawa & Takahashi, 
2006; Hamann et al., 2008b), and we fix the CMB mean tem- 
perature to To — 2.7255 K (Fixsen, 2009). Reionization is as- 
sumed to occur rapidly at a redshift Zre> with optical depth t. 
This model can be characterized by four free cosmological pa- 
rameters: Q.\,h^ ,Q.Ji^ ,6uc,T, defined in Table 1, in addition to 
the parameters describing the initial perturbations. 

3.2. Data 

The CMB data we use consists of the Planck CMB temperature 
likelihood supplemented by the WMAP large-scale polarization 
likelihood (henceforth Planck+WP), as described in Sect. 3.2.1. 
The large scale angular E-mode polarization spectrum is impor- 
tant for constraining reionization. This breaks the degeneracy 
in the temperature data between the primordial power spectrum 
amplitude and the optical depth to reionization. In some cases, 
we also include the Planck lensing likelihood, measurements 
of the CMB at higher resolution from the Atacama Cosmology 
Telescope (ACT) and the South Pole Telescope (SPT), and mea- 
surements of the baryon acoustic oscillation (BAO) scale. This 
follows the same approach as Planck Collaboration XVI (2013). 

3.2.1 . Planck CMB temperature data 

The Planck CMB likelihood is based on a hybrid ap- 
proach, which combines a Gaussian likelihood approximation 
based on temperature pseudo cross-spectra at high multipoles 
(Hamimeche & Lewis, 2008), with a pixel-based temperature 
and polarization likelihood at low multipoles. We summarize the 
likelihood here, and for a detailed description refer the reader to 
Planck Collaboration XV (2013). 
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The small-scale Planck temperature likelihood is based 
on pseudo cross-spectra between pairs of maps at 100, 143, 
and 217 GHz, masked to retain 49%, 31%, and 31% of 
the sky, respectively. This results in angular auto- and cross- 
correlation power spectra covering multipole ranges of 50 < 
{ < 1200 at 100 GHz, 50 < ^ < 2000 at 143 GHz, 
and 500 < £ < 2500 at 217 GHz and for the 143 x 
217 GHz cross-spectrum. In addition to instrumental uncer- 
tainties, mitigated here by using only cross-spectra among dif- 
ferent detectors, small scale foreground and CMB secondary 
anisotropics need to be accounted for. The foreground model 
used in the Planck high-^ likelihood is described in de- 
tail in Planck Collaboration XV (2013); Planck Collaboration 
XVI (2013), and includes contributions to the cross-frequency 
power spectra from unresolved radio point sources, the cos- 
mic infrared background (CIB), and the thermal and kinetic 
Sunyaev-Zeldovich effects. There are eleven adjustable nuisance 

paiameiers. (,Aioo'^143'^217' '^143x217' ^143 '^217 ' M43x217' ^ ' 

). In addition, the calibration parameters for 



tSZ 4kSZ etSZ-CIB 
'143 



the 100 and 217 GHz channels, cioo and Qn, relative to the 
143 GHz channel, and the dominant beam uncertainty eigen- 
mode amplitude Z?j are left free in the analysis, with other beam 
uncertainties marginalized analytically, so the Planck high-^ 
likelihood includes 14 nuisance parameters. 

The low-^" Planck likelihood combines the Planck tempera- 
ture data with the large scale 9-year WMAP polarization data for 
this release. Using the procedure introduced in Page et al. (2007), 
the temperature and polarization likelihood can be separated as- 
suming negligible noise in the temperature map. The temper- 
ature likelihood is based on a Gibbs approach (Eriksen et al., 
2007), mapping out the distribution of the ^ < 50 CMB temper- 
ature multipoles from a foreground-cleaned combination of the 
30 - 353 GHz maps. The polarization likelihood uses a pixel- 
based approach using the WMAP 9-year polarization maps at 
33, 41, and 61 GHz , and includes the temperature-polarization 
cross-correlation (Page et al., 2007). Its angular range is ^ < 23 
for TE, EE, and BB. 



3.2.2. Planck lensing data 

The primary CMB anisotropics are distorted by the gravitational 
potential induced by intervening matter Such lensing, which 
tends to broaden and smooth out the acoustic oscillations, is 
taken into account as a correction to the observed temperature 
power spectrum. The lensing power spectrum can also be recov- 
ered by measuring higher-order correlation functions. 

For some of our analysis we include the Planck lensing 
likelihood, derived in Planck Collaboration XVII (2013), which 
measures the non-Gaussian trispectrum of the CMB, which is 
proportional to the power spectrum of the lensing potential. As 
described in Planck Collaboration XVII (2013), this potential 
is reconstructed using quadratic estimators (Okamoto & Hu, 
2003), and its power spectrum is used to estimate the lensing de- 
flection power spectrum. The spectrum is estimated from the 143 
and 217 GHz maps, using multipoles in the range 40 < ^ < 400. 
The theoretical predictions for the lensing potential power spec- 
trum are calculated at linear order. 



3.2.3. ACT and SPT temperature data 

In some of our analysis we include data from ACT and SPT. 
ACT measures the power spectrum at 148 and 218 GHz, and 
the cross-spectrum (Das et al., 2013), and covers angular scales 



500 < ^ < 10000 at 148 GHz and 1500 < f < 10000 at 218 
GHz. We use these data in the range { > 1000 in combina- 
tion with Planck. SPT measures the power spectrum for angular 
scales 2000 < { < 10000 at 95, 150, and 220 GHz (Reichardt 
et al., 2012). The spectrum at larger scales is also measured at 
150 GHz in Story et al. (2012) but we do not include these in our 
analysis. To model the foregrounds for ACT and SPT we follow 
a similar approach to the likelihood described in Dunkley et al. 
(2013), extending the model used for the Planck high-/" likeli- 
hood. Additional nuisance parameters are included to model the 
Poisson source levels in each experiment, the residual Galactic 
dust, and inter-frequency calibration parameters. They are de- 
scribed in more detail in Planck Collaboration XVI (2013). 

3.2.4. BAO data 

The BAO angular scale acts as a standard ruler and allows us 
to map out the expansion history of the Universe after last scat- 
tering. The BAO scale, extracted from galaxy redshift surveys, 
provides a constraint on the late-time geometry and breaks de- 
generacies with other cosmological parameters. Galaxy surveys 
constrain the ratio Dv(z)/rs, where Dv(z) is the spherically aver- 
aged distance scale to the effective survey redshift z and rs is the 
sound horizon (Mehta et al., 2012). 

In this analysis we consider a combination of the mea- 
surements by the 6dFGRS (Beutler et al. (2011), z = 0.106), 
SDSS-n (Padmanabhan et al. 2012, z = 0.35), and BOSS 
CMASS (Anderson et al. 2012, z = 0.57) surveys, assuming 
no correlation between the three data points. This likelihood is 
described further in Planck Collaboration XVI (2013). 

3.3. Parameter estimation 

Given a model Ai with free parameters x = {xi, ■ ■ ■ , xx} and a 
likelihood function of the data X(data|i), the (posterior) proba- 
bility density !P as a function of the parameters can be expressed 
as 

!P(i|data, M) oc X(data|i) ■ Pix\M), (29) 

where P(x\Ai) represents the data-independent prior probabil- 
ity density. Unless specified otherwise, we choose wide top-hat 
prior distributions for all cosmological parameters. 

We construct the posterior parameter probabilities us- 
ing the Markov Chain Monte Carlo (MCMC) sampler as 
implemented in the CosmoMC (Lewis & Bridle, 2002) or 
MontePython (Audren et al., 2012) packages. In some cases, 
when the calculation of the Bayesian evidence (see below) is 
desired, or when the likelihood function deviates strongly from 
a multivariate Gaussian, we use the nested sampling algorithm 
provided by the MultiNest add-on module (Feroz & Hobson, 
2008; Feroz et al., 2009) instead of the Metropolis-Hastings al- 
gorithm. 

Joint two-dimensional and one-dimensional posterior distri- 
butions are obtained by marginalization. Numerical values and 
constraints on parameters are quoted in terms of the mean and 
68% central Bayesian interval of the respective one-dimensional 
marginalized posterior distribution. 

3.4. Model selection 

Two approaches to model selection are commonly used in statis- 
tics. The first approach examines the logarithm of the likelihood 
ratio, or effective x^. 



Axlff = 2 [ln£„,,(Mi) - lni:„,,(M2)] , 



(30) 
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between models Mi and corrected for the fact that mod- 
els with more parameters provide a better fit due to fitting away 
noise, even when the more complicated model is not correct. 
Various information criteria have been proposed based on this 
idea (Akaike, 1974; Schwarz, 1978), see also Liddle (2007). 
These quantities have the advantage of being independent of 
prior choice and fairly easy to calculate. The second approach 
is a Bayesian approach (Cox, 1946; Jeffreys, 1998; Jaynes & 
Bretthorst, 2003), based on evaluating ratios of the model aver- 
aged likelihood, or Bayesian evidence, defined by 

£, = ^ d}^P{x\Mi)£{AMx). (31) 

Evidence ratios are naturally interpreted as betting odds be- 
tween models.'' Nested sampling algorithms allow rapid numer- 
ical evaluation of &. In this paper we will consider both the ef- 
fective and the Bayesian evidence. 

4. Constraints on slow-roll inflationary models 

In this section we describe constraints on the slow-roll infla- 
tionary paradigm, using Planck+WP data in combination with 
the likelihoods described in Sects. 3.2.2-3.2.4. First we con- 
centrate on characterizing the primordial power spectrum us- 
ing Planck and other data. We start by showing that the em- 
pirical pre-inflationary Harrison-Zeldovich (HZ) spectrum with 
= 1 does not fit the Planck measurements. We further exam- 
ine whether generalizing the cosmological model, for example 
by allowing the number of neutrino species to vary, allowing the 
Helium fraction to vary, or admitting a non-standard reionization 
scenario could reconcile the data with = 1, and conclude that 
Ms 1 is robust. Next we consider evidence for a running of n^ 
and show it constrained to be small, although we find a prefer- 
ence for negative running at modest statistical significance. 

We then investigate the Planck constraints on slow roll in- 
flation, allowing a tilt for the spectral index and the presence of 
tensor modes, and discuss the implications for the simplest stan- 
dard inflationary models. In this section the question is studied 
using the slow roll approximation, but later sections move be- 
yond the slow roll approximation. We show that compared to 
previous experiments, Planck significantly narrows the space of 
possible inflationary models. Finally we comment on the impli- 
cations for inflation of the Planck constraints on possible devia- 
tions from spatial flatness. 

4. 1. Ruling out exact scale invariance 

The simplest ansatz for characterizing the statistical properties 
An of the primordial cosmological perturbations is the so-called 
HZ model proposed by Harrison (1970); Zeldovich (1972); and 
Peebles & Yu (1970). These authors pointed out that a power 
spectrum with exact scale invariance for the Newtonian gravi- 
tation potential fitted the data available at the time, but without 
giving any theoretical justification for this form of the spectrum. 
Under exact scale invariance, which would constitute a sort of 
unexplained new symmetry, the primordial perturbations in the 

' Note that since the average is performed over the entire support 
of the prior probability density, the evidence depends strongly on the 
probability range for the adjustable parameters. Whereas in parameter 
inference, the exact extent of the prior ranges often becomes irrelevant 
as long as they are "wide enough" (i.e., containing the bulk of the high- 
likelihood region in parameter space), the value of the evidence will 
generally depend on precisely how wide the prior range was chosen. 



Newtonian gravitational potential look exactly the same whether 
they are magnified or demagnified. In this simple model, vector 
and tensor perturbations are absent and the spectrum of curva- 
ture perturbations is characterized by a single parameter, the am- 
plitude As. Inflation, on the other hand, generically breaks this 
rescaling symmetry. Although under inflation, scale invariance 
still holds approximately, inflation must end: therefore, as dif- 
ferent scales are being imprinted, the physical conditions must 
also evolve. (However, see Starobinsky (2005) for the canonical 
inflationary model designed to mimic «s = 1, albeit with r + 0). 

Although a detection of a violation of scale invariance would 
not definitively prove that inflation is responsible for the gener- 
ation of the primordial perturbations, ruling out the HZ model 
would confirm the expectation of small deviations from scale in- 
variance, almost always on the red side, which are generic to all 
inflationary models without fine tuning. We examine in detail 
the viability of the HZ model using statistics to compare it to the 
more general model where the spectral index is allowed to vary, 
as motivated by slow-roll inflation. 

When the cosmological model with = 1 is compared with 
a model in which «s is allowed to vary, we find that allowing n^ to 
deviate from one decreases the best-fit effective jt'^ by 27.9 with 
respect to the HZ model. Thus the significance of the finding 
that «s 1 is in excess of 5cr. The parameters and maximum 
likelihood of this comparison are reported in Table 3. 

One might wonder whether n^-\ could be reconciled with 
the data by relaxing some of the assumptions of the underly- 
ing cosmological model. Of particular interest is exploring those 
parameters known to be almost degenerate with the spectral in- 
dex, such as the effective number of neutrino species A^eif and 
the primordial helium fraction }p, which both alter the damp- 
ing tail of the temperature spectrum (Trotta & Hansen, 2004; 
Hou et al., 201 1), somewhat mimicking a spectral tilt. Assuming 
a Harrison-Zeldovich spectrum and allowing A^eff or Ip to float, 
and thus deviate from their standard values, gives almost as good 
a fit to Planck+WP data as the ACDM model with a varying 
spectral index, with A^lg = 2.8 and 2.2, respectively. However, 
as shown in Table 3, the HZ, HZ-nA^eff, and HZ+Yp models re- 
quire significantly higher baryon densities and reionization op- 
tical depths compared to ACDM. In the HZ-nip model, one ob- 
tains a helium fraction of Yp = 0.3194 ± 0.013. This value is in- 
compatible both with direct astrophysical measurements of this 
quantity (Aver et al., 2012) and with standard Big Bang nucle- 
osynthesis (Hamann et al., 2008b). The HZ-nA^eff model, on the 
other hand, would imply the presence of AA^eff ~ 1 new effective 
neutrino species beyond the three known species. When BAO 
measurements are included in the likelihood, A^lg increases to 
39.2 (HZ), 4.6 (HZh-Fp), and 8.0 (HZ-nA^eff), respectively, for the 
three models. The significance of this detection is also discussed 
in Planck Collaboration XVI (2013). 

4.2. Constraining inflationary models using the slow roll 
approximation 

We now consider all inflationary models that can be described 
by the primordial power spectrum parameters consisting of the 
scalar amplitude. As, the spectral index, «s, and the tensor-to- 
scalar ratio r, all defined at the pivot scale We assume that the 
spectral index is independent of the wavenumber k. Negligible 
running of the spectral index is expected if the slow roll con- 
dition is satisfied and higher order corrections in the slow roll 
approximations can be neglected. In the next subsection we re- 
lax this assumption. 
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HZ 


HZ + ¥p 


HZ + Weff 


ACDM 


100 Smc 

T 

^(10'%) 
Na 


2296 ± 24 
1088+ 13 
1.04292 + 0.00054 

0.125!«;5!|5 

T 1 TQ+0.032 

-'•^-'-'-0.028 


2296 + 23 
1158 + 20 
1.04439 + 0.00063 
0.109!!J;!!|3 

■3 1 Q7+0.027 
-^ ■-^ -0.028 

0.3194 + 0.013 


2285 + 23 
1298 + 43 
1.04052 + 0.00067 
0-105t°|^ 

T 140+0.027 

3.98 + 0.19 


2205 + 28 
1199 + 27 
1.04131 + 0.00063 

3 089+""^'* 
0.9603 + 0.0073 


-2Aln(£„,,) 


27.9 


2.2 


2.8 






Table 3. Constraints on cosmological parameters and best-fit -2Aln(X) with respect to the standard ACDM model, using 
Planck+WP data, testing the significance of the deviation from the HZ model. 



Sampling the power spectrum parameters A^, rig, and r is 
not the only method for constraining slow roll inflation. Another 
possibility is to sample the HFF in the analytic expressions for 
the scalar and tensor power spectra (Stewart & Lyth, 1993; Gong 
& Stewart, 2001; Leach et al., 2002). In the Appendix, we per- 
form a comparison of slow-roll inflationary predictions by sam- 
pling the HFF with Planck data, and show that the results ob- 
tained in this way agree with those derived by sampling the 
power spectrum parameters. This confirms similar studies with 
previous data (Hamann et al., 2008c; FinelU et al., 2010). 

The spectral index estimated from Planck+'WP data is 



n, = 0.9603 + 0.0073. 



(32) 



This tight bound on is crucial for constraining inflation. The 
Planck constraint on r depends slightly on the pivot scale; we 
adopt kt = 0.002 Mpc"' to quote our result, with ro ooa < 0.12 
at 95% CL. This bound improves on the most recent results, 
including the WMAP 9-year constraint of r < 0.38 (Hinshaw 
et al., 2012a), the WMAP7 + ACT Hmit of r < 0.28 (Sievers 
et al., 2013), and the WMAP7 + SPT Hmit of r < 0.18 (Story 
et al., 2012). The new bound from Planck is consistent with 
the limit from temperature anisotropics alone (Knox & Turner, 
1994). When a possible tensor component is included, the spec- 
tral index from Planck+WP is not significantly changed, with 
n, = 0.9624 ± 0.0075. 

The Planck constraint on r corresponds to an upper bound 
on the energy scale of inflation 



rM^i = (1.94 X 10'" GeV> 



0.12 ' 



(33) 



at 95% CL. This is equivalent to an upper bound on the Hubble 
parameter during inflation of Ht/Mpi < 3.7 x 10"^. In terms of 
slow-roll parameters, Planck+WP constraints imply ey < 0.008 
at 95% CL, and 7]v = -0.010!^;™^^. 

The Planck results on «s and r are robust to the addition 
of external data sets (see Table 4). When the high-/" CMB 
ACTh-SPT data are added, we obtain n, = 0.9600 ± 0.0072 and 
'"0.002 < 0.11 at 95% CL. Including the Planck lensing likeli- 
hood gives Hs = 0.9653 ± 0.0069 and ro.002 < 0.13, and adding 
BAO data gives n, = 0.9643 ± 0.0059 and ro.002 < 0.12. These 
bounds are robust to the small changes in the polarization likeli- 
hood at low multipoles. To test this robustness, instead of using 
the WMAP polarization likelihood, we impose a Gaussian prior 
T = 0.07 + 0.013 to take into account small shifts due to un- 
certainties in residual foreground contamination or instrument 
systematics in the evaluation of t, as performed in Appendix B 
of Planck Collaboration XVI (201 3). We find at most a reduction 
of 8% for the upper bound on r. 



It is useful to plot the inflationary potentials in the n^-r plane 
using the first two slow-roll parameters evaluated at the pivot 
scale k^ = 0.002 Mpc"' (Dodelson et al., 1997). Given our ig- 
norance of the details of the epoch of entropy generation, we 
assume that the number of e-folds to the end of inflation lies 
in the interval [50,60]. This uncertainty is plotted for those po- 
tentials predicting an exit from inflation without changing the 
potential. 

Fig. 1 shows the Planck constraints in the - r plane and 
indicates the predictions of a number of representative inflation- 
ary potentials. The sensitivity of Planck data to high multipoles 
removes the degeneracy between n^ and r found using WMAP 
data. Planck data favour models with a concave potential. As 
shown in Fig. 1, most of the joint 95% allowed region lies be- 
low the convex potential limit, and concave models with a red 
tilt in the range [0.945-0.98] are allowed by Planck at 95% CL. 
In the following we consider the status of several illustrative and 
commonly discussed inflationary potentials in light of the Planck 
observations. 

Power law potential and chaotic inflation 

The simplest class of inflationary models is characterized by a 
single monomial potential of the form 



Pi M, 



(34) 



This class of potentials includes the simplest chaotic models, in 
which inflation starts from large values for the inflaton, > Mpi. 
Inflation ends by violation of the slow-roll regime, and we as- 
sume this occurs at ey = 1. According to Eqs. 5, 6, and 15, 
this class of potentials predicts to lowest order in slow-roll pa- 
rameters «s - 1 -n(n + 2)Mp,/(^J, r ^ 8n'^M^J4>l, (pi » 

nMpj(4A^, + n)/2. The A<p'^ model lies well outside of the joint 
99.7% CL region in the n^ - r plane. This result confirms pre- 
vious findings from e.g., Hinshaw et al. (2012a) in which this 
model is well outside the 95% CL for the WMAP 9-year data 
and is further excluded by CMB data at smaller scales. 

The model with a quadratic potential, n - 2 (Linde, 1983), 
often considered the simplest example for inflation, now lies 
outside the joint 95% CL for the Planck+WP +h\gh-£ data for 
Nt < 60 e-folds, as shown in Fig. 1. 

A linear potential with n - I (McAllister et al., 2010), mo- 
tivated by axion monodromy, has rjv - Q and lies within the 
95% CL region. Inflation with n - Ijli (Silverstein & Westphal, 
2008), however, also motivated by axion monodromy, now lies 
on the boundary of the joint 95% CL region. More permissive 
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Model 


Parameter 


Flanck+WP 


P/ancfc+WP+lensing 


Planck + WP+high-^" 


Planck+WP+BAO 


ACDM + tensor 




0.9624 ± 0.0075 
< 0.12 


0.9653 ± 0.0069 
< 0.13 


0.9600 ±0.0071 
< 0.11 


0.9643 + 0.0059 
< 0.12 


-lAlnXn,^ 











-0.31 



Table 4. Constraints on the primordial perturbation parameters in the ACDM+r model from Planck combined with other data sets. 
The constraints are given at the pivot scale k^ - 0.002 Mpc"'. 
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Fig. 1. Marginalized joint 68% and 95% CL regions for and ro.oo2 from Planck in combination with other data sets compared to 
the theoretical predictions of selected inflationary models. 



reheating priors allowing Nt < 50 could reconcile this model 
with the Planck data. 

Exponential potential and power law inflation 
Inflation with an exponential potential 

y(<^) = A^expJ-^^j (35) 

is called power law inflation (Lucchin & Matarrese, 1985), 
because the exact solution for the scale factor is given by 
a(t) OC f^/'*'. This model is incomplete, since inflation would 
not end without an additional mechanism to stop it. Assuming 
such a mechanism exists and leaves predictions for cosmo- 
logical perturbations unmodified, this class of models predicts 
r - -8(ns - 1) and is now outside the joint 99.7% CL contour 

Inverse power law potential 

Intermediate models (Barrow, 1990; Muslimov, 1990) with in- 
verse power law potentials 

V(^) = A^(^f (36) 



lead to inflation with a(t) oc exp(Af^), with A > and < / < 1, 
where / - 4/(4 + /3) and > 0. In intermediate inflation there 
is no natural end to inflation, but if the exit mechanism leaves 
the inflationary predictions on cosmological perturbations un- 
modified, this class of models predicts r » -8/3(«s - - 2) 
(Barrow & Liddle, 1993). It is disfavoured, being outside the 
joint 95% CL contour for any j6. 

Hill-top models 

In another interesting class of potentials, the inflaton rolls away 
from an unstable equilibrium as in the first new inflationary mod- 
els (Albrecht & Steinhardt, 1982; Linde, 1982). We consider 

y(^).A4jl_^ + ...j, (37) 

where the elUpsis indicates higher order terms negligible during 
inflation, but needed to ensure the positiveness of the potential 
later on. An exponent of /? = 2 is allowed only as a large field 
inflationary model and predicts n^ - 1 ~ —AM^J^^ + 3r/2i and 

r X! 320jMpj/yu'*. This potential leads to predictions in agree- 
ment with Planck+WP+BAO joint 95% CL contours for super- 
Planckian values of jj, i.e., > 9 Mpi. 

Models with p >3 predict n^ - I ~ -(2/N)(p - l)/(p - 2) 
when r ~ 0. The hill-top potential with p - 3 lies outside the 
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joint 95% CL region for Planck+WP+BAO data; the case with 
p - 4 is also in tension with P/fl«cfe+WP+BAO, but allowed 
within the joint 95% CL region for Nt > 50. For larger values of 
r these models provide a better fit to the Planck+'WP+BAO data. 

A simple symmetry breaking potential 

The symmetry-breaking potential (Olive, 1990) 



A' I 



1^' 



(38) 



can be considered as a self-consistent completion of the hill-top 
model with p - 2 (although it has a different limiting large- 
field branch for non-zero r). This potential leads to predictions 
in agreement with Planck + WP + B AO joint 95% CL contours 
for super-Planckian value of jj, i.e. fj.> 13 Mpi. 

Natural inflation 

Another interesting class of potentials is natural inflation 
(Freese et al., 1990; Adams et al., 1993), initially motivated by 
its origin in symmetry-breaking in an attempt to naturally give 
rise to the extremely flat potentials required for inflationary cos- 
mology. In natural inflation the effective one-dimensional po- 
tential takes the form 



V{<p) = A^ 



1 + cos 



(39) 



where / is a scale which determines the slope of the potential. 
Depending on the value of /, the model falls into the large field 
(/ > 1.5 Mpi) or small field (/ < 1.5 Mpi) classification scheme. 
Therefore, n, ^ I - M^^/f holds for small / and ~ 1 - 2/N, 

r X ^/N holds for large /, approximating the m^<p^ potential in 
the latter case (with A^, ^ (2f/M^{)ln[sm((pJf)/sm((pJf)]). 
This model agrees with Planck+WP data for / > 5 Mp]. 



predicts ns - I ~ -(1 + 3a;h/2)/A^t and r ^ Sffh/A^*. For aj, 1 
and A^, ^ 50, n, =^ 0.98 is disfavoured by Planck+WP+BAO 
data at more than 95% CL. However, more permissive reheating 
priors allowing Nt < 50 or a non-negligible ah give models that 
are consistent with the Planck data. 



inflation 

Inflationary models can also be accommodated within extended 
theories of gravity. These theories can be analysed either in the 
original (Jordan) frame or in the conformally-related Einstein 
frame with a Klein-Gordon scalar field. Due to the invariance of 
curvature and tensor perturbation power spectra with respect to 
this conformal transformation, we can use the same methodol- 
ogy described earlier 

The first inflationary model proposed was of this type 
and based on higher order gravitational terms in the action 
(Starobinsky, 1980) 



S=/..,V=J^(«.^), (43) 

with the motivation to include semi-classical quantum effects. 
The predictions for inflation were first studied in Mukhanov 
& Chibisov (1981); Starobinsky (1983) and can be summarized 
as: n,-l x -8(4A^,-h9)/(4A?.h-3) and r « 192/(4A^.+3)2. Since r 
is suppressed by another l/N, with respect to the scalar tilt, this 
model predicts a tiny amount of gravitational waves. This model 
predicts ns - 0.963 for - 55 and is fully consistent with the 
Planck constraints. 



Non-minimally coupled inflaton 

A non-minimal coupling of the inflaton to gravity with the action 



Hybrid inflation 

In hybrid inflationary models a second field, coupled to the 
inflaton, undergoes symmetry breaking. The simplest example 
of this class is 



y(0,;r) = A^(l-^ 



+ UW + |-0Y . (40) 



For most of their parameter space, these models can be consid- 
ered effectively as single field models for the inflaton <p. The 
second field x is close to the origin during the slow -roll regime 
for 0, and inflation ends either by breakdown of slow roll for 
the inflaton at ^ (df//d0)V(A'* + U{^)f ~ 1 or by the 
waterfall transition of;^^. The simplest models with 



U{<P) 



(41) 



are disfavoured for most of the parameter space (Cortes & 
Liddle, 2009). Models with m^(jp-l2 ~ A* are disfavoured due to 
a high tensor-to-scalar ratio, and models with U ((p) <K A'* predict 
a spectral index «s > 1. also disfavoured by the Planck data. 

We discuss hybrid inflationary models predicting n^ < \ sep- 
arately. As an example, the spontaneously broken SUSY model 
(Dvah et al., 1994) 



U{(p) = ffhA'* In 



(42) 



(44) 

leads to several interesting consequences such as a lowering of 
the tensor-to-scalar ratio. 

The case of a massless self-interacting inflaton (0o - 0) 
agrees with the Planck+WP data for ^ Q. Within the range 
50 < < 60, this model is within the Planck+WP joint 
95% CL region for ^ > 0.0019, improving previous bounds 
(Okada et al., 2010). 

The amplitude of scalar perturbations is proportional to A/^^ 
for ^ » 1, and therefore the problem of tiny values for the in- 
flaton self-coupling A can be alleviated (Salopek et al., 1989; 
Fakir & Unruh, 1990). The regime <po <s Mpi is allowed and <p 
could be the Standard Model Higgs as proposed in Bezrukov & 
Shaposhnikov (2009). The Higgs case with ^ s> 1 has the same 
predictions as the R^ model in terms of «s and r as a function 
of Nf The reheating mechanism in the Higgs case can be more 
efficient than in R^ case and therefore predicts a slightly larger 
ns- This model is fully consistent with Planck constraints. 

The case with ^ < and \^\4>l/M^^ ~ 1 was also recently 
emphasized in Linde et al. (2011). With the symmetry breaking 
potential in Eq. 44, the large field case with inflaton </> > 0o is 
disfavoured by Planck data, whereas the small field case (p < (po 
is in agreement with the data. 
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4.3. Running spectral index 

We have shown that the single parameter Harrison-Zeldovich 
spectrum does not fit the data, and that at least the first two terms. 
As and n^, in the expansion of the primordial power spectrum in 
powers of ln{k) given in Eq. 10 are needed. Here we consider 
whether the data require the next term known as the running of 
tlie spectral index (Kosowsky & Turner, 1995), defined as the 
derivative of the spectral index with respect to Infc, dnj.t/dlnfc 
for scalar or tensor fluctuations. In slow-roll single-field infla- 
tion, the running is second order in the Hubble slow-roll param- 
eters, for scalar and for tensor perturbations (see Eqs. 17 and 
18, respectively) (Leach et al., 2002), and thus is typically sup- 
pressed with respect to, e.g., - 1 and n^, which are first order 
If the slow-roll approximation holds and the inflaton has reached 
its attractor solution, d«s/dlnfc and dnt/dlnA; are related to the 
potential slow -roll parameters as in Eqs. 17 and 18. Given the 
tight constraints on the first two slow-roll parameters ey and rfy 
(e\ and 62) from the present data, typical values of the running to 
which Plancic is sensitive (Pahud et al., 2007) would generically 
be dominated by the contribution from the third derivative of the 
potential, encoded in \^y\ (or 63). 

While it is easy to see that the running is invariant under a 
change in pivot scale, the same does not hold true for the spectral 
index and the amplitude of the primordial power spectrum. It is 
convenient to choose ^» such that dWs/dlnA: and n^ are decorre- 
lated (Cortes et al., 2007). This approach minimizes the inferred 
variance of n^ and facilitates comparison with constraints on n^ 
in the power law models. Note however that the decorrelation 
pivot scale Icf^'^ depends both on the model and data set consid- 
ered in the analysis. 

We consider a model parameterizing the power spectrum us- 
ing As(A:,) , ns(A:,), and dns/dlnfc, where A:, - 0.05 Mpc"'. The 
joint constraints on n^ and dn^/dln Ic at the decorrelation scale 
of /tf = 0.04 Mpc"' are shown in Fig. 4. The PlancIc+WP con- 
straints on the running do not change significantly when com- 
plementary data sets such as Plancic lensing, CMB high-/", and 
BAO data are included. We find 

dnjdlnlc = -0.013 ± 0.009 (68% CL, Planclc+WP) , (45) 

which is negative at the 1.5cr level. This reduces the the uncer- 
tainty compared to previous CMB results. Error bars are reduced 
by 60% compared to the WMAP 9-year results (Hinshaw et al., 
2012a), and by 20-30% compared to WMAP supplemented by 
SPT and ACT data (Hou et al., 2012; Sievers et al., 2013). Plancic 
finds a smaller scalar running than SPTH-WMAP7 (Hou et al., 
2012), and larger then ACTH-WMAP7 (Sievers et al., 2013). The 
best-fit likelihood improves by only ^xlg ~ 1-5(3 when high-/" 
data are included) with respect to the minimal case in which is 
scale independent, indicating that the deviation from scale inde- 
pendence is not very significant. The constraint for the spectral 
index in this case is 0.9630 + 0.0065 at 68% CL at the decor- 
relation pivot scale A:, = 0.038 Mpc"'. This result implies that 
the third derivative of the potential is small, i.e., \^y\ ~ 0.007, 
but compatible with zero at 95% CL, for inflation at low energy 
(i.e., with ey ~ 0). 

We also test the possibility that the running depends on the 
wavelength so that d^Ws/dlnA:^ is nonzero. With Planclc+WP 
data, we find d-ns/dlnA:^ = 0.020+°°[^. This resuh is stable 
with respect to the addition of complementary data sets. When 
d^«s/dlnA;^ is allowed in the fit, we find a value for the running 
dws/d In k consistent with zero. 

Finally we allow a non-zero primordial gravitational wave 
spectrum together with the running. The tensor spectral in- 
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Fig. 2. Marginalized joint 68% and 95% CL for (dnjdln k,n,) 
using Planck+WP+BAO, either marginalizing over r or fixing 
r = at A:* = 0.038 Mpc"' . The purple strip shows the prediction 
for single monomial chaotic inflationary models with 50 < A^* < 
60 for comparison. 



dex and its running are set by the slow-roll consistency re- 
lations to second order, with «( = -'"(2 - ?"/8 - ns)/8 and 
dnjdlnk - r(r/8 + Ws ~ l)/8. Planck measures the running to 
be dn^/dlnk - -0.016 ± 0.010 when tensors are included (see 
Table 5 and Fig. 4). The constraints on the tensor-to-scalar ra- 
tio are relaxed compared to the case with no running, due to an 
anti-correlation between r and dn^/dlnk, as shown in Fig. 4 for 
Planck+WP+BAO. 

Varying both tensors and running, Planck+WP implications 
for slow-roll parameters are ey < 0.015 at 95% CL, rjv = 
-0.014!^;°j5^ \^l\ = 0.009 ± 0.006. 

In summary, the Planck data prefer a negative running for 
the scalar spectral index of order dn^/dlnk ^ -0.015 at the 
1.5cr significance level, alone and in combination with other 
astrophysical data sets. Weak statistical evidence for non-zero 
negative values of dn^/dlnk has been claimed in several previ- 
ous investigations with the WMAP data and smaller scale CMB 
data (e.g., Spergel et al., 2003; Peiris et al., 2003; Dunkley et al., 
2011; Hinshaw et al., 2012a; Hou et al., 2012). 

If primordial, negative values for dn^/dlnk of order 10"^, 
would be interesting for the physics of inflation. The running of 
the scalar spectral index is a key prediction for inflationary mod- 
els. It is strictly zero for power law inflation, whose fit to Planck 
was shown to be quite poor in the previous section. Chaotic 
monomial models with V((p) oc 0" predict d«s/dlnA: x -8(« + 
2) / (4A^ + n)^ a: (hs - 1 )^, and the same order of mag nitude(10"^) 
is quite typical for many slow-roll inflationary models, such as 
natural inflation or hilltop inflation, to name a few. It was pointed 
out that a large negative running of d«s/dlnA; < -10"^ would 
make it difficult to support the A^, ^ 50 e-foldings required from 
inflation (Easther & Peiris, 2006), but this holds only without 
nonzero derivatives higher than the third order in the inflation- 
ary potential. Designing inflationary models that predict a neg- 
ative running of 0(10^^) with an acceptable and number of 
e-folds is not impossible, as the case with modulated oscilla- 
tions in the inflationary potential demonstrates (Kobayashi & 
Takahashi, 2011). This occurs, for instance, in the axion mon- 
odromy model when the instanton contribution is taken into ac- 



Planck Collaboration: Constraints on inflation 



13 



Model 


Parameter 


Planck+WP 


P/ancfc+WP+lensing 


Planck+WP+high-£ 


/'/ancvt+WP+BAO 






0.9561 ± 0.0080 


0.9615 ±0.0072 


0.9548 ± 0.0073 


0.9596 ± 0.0063 


ACDM + dnjdlnk 


dils/dlnfc 


-0.0134 ±0.0090 


-0.0094 ± 0.0085 


-0.0149 ±0.0085 


-0.0130 ±0.0090 




-2Aln T 


-1.50 


-0.77 


-2.95 


-1.45 


ACUM + an,/Qln/c 


dns/dlnx 
d^iijdlnfc^ 


n r\r\ i +0.016 
"•"■^"-0.015 


n Q'^7^+0 077 

n nn^+0.015 
0-006_ooi4 
019+0018 


0.9476_oo88 
n nni +0.013 
0.001_ooi4 
022+" '"'' 

U.UiZ._oo[3 


n 0^^^8+0.068 

0.9568 

n nnn+0.0l6 

0.000_ooi3 
017+° ""' 

U.Ul /_ooi4 




-9A In r 




-2.14 




-? 4.0 






0.9583 ±0.0081 


0.9633 ± 0.0072 


0.9570 ± 0.0075 


0.9607 ± 0.0063 




r 


< 0.25 


<0.26 


< 0.23 


< 0.25 


ACDM + r + dnjdlnk 


d7is/dlnfc 


0.021 ±0.012 


0.017 ±0.012 


-0 022+'' °" 


-0 021+° °'- 

"■"^^+0.010 




-2Alni;,^ax 


-1.53 


-0.26 


-3.25 


-1.5 



Table 5. Constraints on the primordial perturbation parameters for ACDM-ndns/dlnfc, ACDM-ndns/d In A;H-r and 
ACDMH-dns/dlnA:H-d^ns/dlnfc^ models from Planck combined with other data sets. Constraints on the spectral index and its de- 
pendence on the wavelength are given at the pivot scale of k^ - 0.05 Mpc 
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Fig. 3. Marginalized joint 68% and 95% CL regions for 
(d^nj din k^ , dnjdlnk) using Planck+WP+BAO. 



Fig. 4. Marginalized joint 68% and 95% CL regions for (r ,n^), 
using P/ancfe-i-WP-i-BAO with and without a running spectral in- 
dex. 



count (McAllister et al., 2010) giving the potential 

y(0) A'^cos(^) . (46) 



4.4. Open inflation 

Most models of inflation predict a nearly flat spatial geome- 
try with small deviations from perfect spatial flatness of order 
\Q.k\ ~ 10"^. Curvature fluctuations may be regarded as local 
fluctuations in the spatial curvature, and even in models of infla- 
tion where the perturbations are calculated about a spatially flat 
background, the spatial curvature on the largest scales accessible 
to observation now are subject to fluctuations from perfect spa- 
tial flatness (i.e., Q.^ = 0). This prediction for this fluctuation is 
calculated by simply extrapolating the power law spectrum to the 
largest scale accessible today, so that Q.^ as probed by the CMB 
roughly represents the local curvature fluctuation averaged over 
our (causal) horizon volume. Although it has sometimes been 
claimed that spatial flatness is a firm prediction of inflation, it 



was realized early on that spatial flatness is not an inexorable 
consequence of inflation, and large amounts of spatial curvature 
(i.e., large compared to the above prediction) can be introduced 
in a precise way while retaining all the advantages of inflation 
(Gott, 1982; Gott & Statler, 1984) through bubble nucleation by 
false vacuum decay (Coleman & De Luccia, 1980). This pro- 
posal gained credence when it was shown how to calculate the 
perturbations in this model around and beyond the curvature 
scale (Bucher et al., 1995; Bucher & Turok, 1995; Yamamoto 
et al., 1995; Tanaka & Sasaki, 1994). See also (Ratra & Peebles, 
1995, 1994; Lyth &. Stewart, 1990). For more refined later cal- 
culations see for example Garriga et al. (1998, 1999); Gratton & 
Turok (1999) and references therein. For predictions of the ten- 
sor perturbations see for example Bucher & Cohn (1997); Sasaki 
et al. (1997); Hertog & Turok (2000). 

An interesting proposal using singular instantons and not 
requiring a false vaccum may be found in Hawking & Turok 
(1998), and for calculations of the resulting perturbation spectra 
see (Hertog & Turok, 2000; Gratton et al., 2000). Models of this 
sort have been studied more recently in the context of the string 
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Fig. 5. Joint posterior for tensors and running of the scalar 
spectral index marginalizing over other parameters. The con- 
tours are set at 68% and 95%. The red contours apply for 
f/flnc^+ WP+high-^ data. The colour of the scattered points in- 
dicates the distribution of Q.^ in this case, he blue contours apply 
when BAO data is also included. (Q/f is then found to be well- 
constrained close to zero.) Solid black hnes guide the eye to the 
no-tensors, no-running solution. 



landscape (see for example Vilenkin (2007) for a nice review). 
Although some proposals for universes with positive curvature 
within the framework of inflation have been put forth (Gratton 
et al., 2002), it is much harder to obtain a closed universe with 
a spatial geometry of positive spatial curvature (i.e., Q.^ < 0) 
(Linde, 2003). 

Theoretically, it is of interest to measure Q.k to an accuracy 
of approximately 10""* or slightly better to test the extrapolation 
of simple flat inflation for this observable. A statistically signif- 
icant positive value would suggest that open inflation, perhaps 
in the context of the landscape, was at play. A statistically sig- 
nificant negative value could pose difficulties for the inflationary 
paradigm. For a recent discussion of these questions, see for ex- 
ample Freivogel et al. (2006), Kleban & Schillo (2012) and Guth 
& Nomura (2012). 

In order to see how much spatial curvature is allowed, 
we consider a rather general model including the param- 
eters r, Hs, and dn^/dlnk as well as Q.k- We find that 

= -0.058+JJ{J^^ with Planck+WP, and = -0.0004 ± 
0.00036 with P/fl«cA:H-WPH-BAO. More details can be found by 
consulting the parameter tables available online. Figure 5 shows 
r and d«s/dlnfc for this family of models. We conclude that any 
possible spatial curvature is small in magnitude even within this 
general model and that the spatial curvature scale is constrained 
to lie far beyond the horizon today. Open models predict a ten- 
sor spectrum enhanced at small wavenumber k<l, where ^ = 1 
corresponds to the curvature scale, but our constraint on Q.^ and 
cosmic variance imply that this aspect is likely unobservable. 



* Available at: http : / / www . sciops .esa. int/ index . php 
?pro ject=planck&page=Planck_LegacY_Ar chive 



4.5. Relaxing tlie assumption of tlie late-time 
cosmological concordance model 

The joint constraints on and r shown in Fig. 1 are one of the 
central results of this paper. However, they are derived assum- 
ing the standard ACDM cosmology at late times (i.e., z < lO'*). 
It is therefore natural to ask how robust our conclusions are to 
changes of the late time cosmological model. We discuss two 
classes of models: firstly, changes to the ACDM energy content, 
and secondly, a more general reionization model. These exten- 
sions can lead to degeneracies of the additional parameters with 
«s or r. 

4.5.1 . Extensions to the energy content 

We consider the ACDM-nrH-A^eff, ACDMH-rH-Fp, ACDMH-rH-^ niy 
and ACDM+r+w extensions of the standard model. This se- 
lection is motivated by the impact on the CMB damping tail 
of the first two and the effect on the Sachs-Wolfe plateau at 
low multipoles for the latter two. The resulting contours are 
shown in Fig. 6. While the lower limit on is stable under 
all extensions considered here, the models that affect the high- 
{ part of the spectrum permit significantly bluer spectral tilts, 
and accordingly also leads to a weaker bound on the tensor-to- 
scalar ratio (ACDM-HrH-A^eff gives n^ - 0.9764 ± 0.0106 and 
roooi < 0.15 at 95% CL, while for ACDM-nrH-yp we obtain 
n, = 0.9810 ± 0.01 1 1 and ro.oo2 < 0.18 at 95% CL). The models 
affecting the large-scale part of the power spectrum, on the other 
hand, do not lead to a notable degradation of constraints on ei- 
ther «s or r (n, = 0.9648 ± 0.0061 and rooo2 < 0.13 at 95% CL 
for ACDM+r+Y, niy, and n, = 0.9601 + 0.0070 and ro.oo2 < 0.11 
at 95% CL for ACDM+r+w). 

4.5.2. General reionization scenario 

In the standard rapid reionization scenario typically used in 
CMB analysis, the Universe is assumed to be completely trans- 
parent after recombination, but the ionization fraction increases 
from zero to one over a duration Az ~ 1 at a certain redshift 
Zreiom which is the only unknown parameter of the reionization 
model. This model is obviously simplistic, but for CMB analy- 
sis it works quite well because the CMB has little sensitivity to 
the details of how the ionization fractrion changes from to 1 . 
In this section we study to what extent allowing more general 
reionization scenarios may alter some of the conclusions con- 
cerning the constraints on and r as well as r. As discussed in 
Mortonson & Hu (2008a) and Mortonson & Hu (2008b), CMB 
anisotropics constrain reionization almost entirely by using the 
shape of the large-scale EE power spectrum, and the power is 
redistributed from larger to smaller scales for reionization pro- 
cesses which take place during a non-negligible redshift interval, 
since they start at an earlier epoch. 

We use the method developed by Mortonson & Hu (2008a) 
to describe and constrain the reionization history. A complete 
principal component basis serves to describe the effect of reion- 
ization on the large-scale E-mode polarization power spectrum. 
Following Mortonson & Hu (2008a) we bin the ionization his- 
tory Xe(z,) using 95 equal width bins with Az - 0.25 ranging 
from Zmin = 6 to Ztnax = 30. For the redshifts z < Zmin we 
assume values for Xe which take into account first (and possi- 
bly second) helium ionization and complete hydrogen ioniza- 
tion (xg = 1.16 for z < 3 and - 1.08 for z < 6). For z > 30 
we fix Xe = 2 X 10""* as the value of Xg expected before reion- 
ization (and after primordial recombination). Any reionization 
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Fig. 6. Marginalized joint 68% and 95% CL regions for 
Planck+WP+BAO data for ACDM+r+A^eff and ACDM+r+Fp 
(top); and ACDM+r+2 my and ACDM+r+w (bottom). Shown 
for comparison are the corresponding contours for the ACDM+r 
model. 

history can be parametrized as a free function of redshift by de- 
composing the ionization fraction as Xe(z) - xl(z) + 2^ m^^S ^^(z), 
where the principal components, 5^(z), are the eigenfunctions 
of the Fisher matrix computed by taking the derivatives of the 
EE polarization power spectrum with respect to x^(z). Following 
Mortonson & Hu (2008a) we consider here the first five eigen- 
functions, S ^(z) with fi - I, ...,5, which will be varied with the 
other cosmological parameters.. 

In Fig. 7 we plot the 68% and 95% CL regions for «s and 
r. The constraint on the tensor-to-scalar ratio is not signifi- 
cantly affected by this additional marginalization, increasing to 
'"0.002 < 0.13 at 95% CL. The scalar spectral index is increased 
to = 0.9650 + 0.0080, compared with n, = 0.9603 ± 0.0073 
obtained with the rapid reionization scheme. This is the same 
trend as noted in Pandolfi et al. (2010) using WMAP data, but 
the effect is less significant due to the improved measurement 
of the spectrum by Planck. The larger freedom in the reioniza- 
tion history increases the posterior on the derived optical depth, 
which is still partially degenerate with the scalar spectral index 
of primordial perturbations. The - \ model is still excluded 
at high significance; however, we find Ax\fi - 12.5 compared to 
the ACDM model. 
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Fig. 7. Marginahzed joint 68% and 95% CL regions for 
Planck+WP data for ACDM-i-r for instantaneous and general 
reionization. 

5. Inflationary model comparison 

In Sect. 4.2 several representative families of parameterized 
models for the inflationary potential were analysed within the 
slow-roll approximation in the neighbourhood of the pivot scale 
kt. Approximate constraints were applied to reject models for 
which there is no plausible scenario for entropy generation. In 
this section we revisit some of the parametric models defined in 
Sect. 4.2. Here, however, the modes for the first order perturba- 
tions, as described in Sect. 2.2, are integrated numerically. Thus 
there is no slow-roll approximation, and the issue of the exis- 
tence of a plausible scenario for entropy generation is examined 
more carefully. We perform a statistical model comparison be- 
tween the competing parameterized potentials, both within the 
framework of Bayesian model comparison and in terms of the 
relative likelihoods of the best-fit models from each parameter- 
ized family. 

As noted in Sect. 2.3, there is considerable uncertainty sur- 
rounding what occurs during the epoch of entropy generation, 
particularly with respect to the energy scale at which entropy 
generation ends, and the average equation of state between that 
epoch and the end of inflation. For this reason, we explore 
a number of scenarios for this intermediate era during which 
entropy generation takes place. The models compared in this 
section include inflation with power law potentials, defined in 
Eq. 34, with several representative values for n; natural inflation, 
in Eq. 39; and hilltop inflation, with p - A and A - 4A'^/f/ in 
Eq. 37. 

The free parameters in these inflationary potentials may vary 
over several orders of magnitude corresponding to unknown 
scales in high energy particle physics. Consequently a logarith- 
mic prior is a sensible choice for these parameters. However, 
there is no theoretical guidance on how to truncate these model 
priors. We therefore adopt broad priors initially, and then trun- 
cate them as follows in order to compare all the models on an 
equal footing. 

The strongest constraint on the inflationary parameter space 
comes from the amplitude of the primordial power spectrum. 
This is a free parameter in most models, and successful structure 
formation in a universe dominated by cold dark matter has long 
been known to require primordial fluctuations with fi ^ 10"^, 
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n = 2/3 


-13<logioW<-7 


n = 1 


-13<logioW<-7 


11 = 2 


-13.5 < log,o(m2) < -8 


« = 4 


-16<log,oW<-10 


1 'I ctLUI tXL 


— S <^ 1n(T (A") <: 




< log,o(/) < 2.5 


Hilltop 


-8<logio(A)<-l 


p = 4 


-17<log,oW<-10 


ACDM 


0.9 < n, < 1.02 




3.0 < logio[lO"'A,] < 3.2 



Matching 


Prior 


A'. 


20 < W < 90 



Table 6. Model priors. Dimensionful quantities are expressed in 
units with Mpi set to unity. The matching parameter A^*, com- 
puted at - 0.05 Mpc"', allows us to marginalize over the 
uncertainty in connecting the inflationary era to astrophysical 
scales. Additional cuts are made in these parameter ranges to se- 
lect out physically relevant subspaces which are compatible with 
basic structure-formation requirements and different entropy- 
generation scenarios, as described in the text, resulting in non- 
rectangular model priors in some cases. We also marginalize 
over the concordance cosmological parameters and foreground 
parameters of the Planck likelihood, as described in the text. 



or A,s K. 10">" (see e.g., Zeldovich, 1972; Linde, 1990). We can 
therefore immediately reject models for which As is far from 
this value, so regions of parameter space which do not yield 
10"" < As < 10"^ are a priori excluded. This range is gen- 
erous relative to estimates of As prior to Planck (e.g., Komatsu 
et al., 2011), but the results do not depend strongly on the 
range chosen. This effectively truncates the logarithmic priors 
on the model parameters, leaving a parameter subspace compat- 
ible with basic structure formation requirements. 

For the single parameter models this requirement defines the 
range of A in Eq. 34. However, for generic multi-parameter mod- 
els, an As-based cut may select a nontrival region of parameter 
space, as happens for the two cases considered here (see e.g.. 
Fig. 1 in Easther & Peiris (2012)). Without the As-based cut in 
the prior, the parameter volume for both natural and hilltop infla- 
tion would be rectangular, and the corresponding Bayesian evi- 
dence values computed for these models would be lowered. 

As discussed in Sect. 2.3, specifying an inflationary poten- 
tial does not enable us to predict the late-time CMB angular 
power spectra. The subsequent expansion history and details of 
the epoch of entropy generation are required to relate the value 
of the inflaton field at horizon crossing to comoving wavenum- 
bers in today's Universe, via Eq. 24. Physically, the fundamental 
parameter that sets the observable perturbation spectrum is the 
value of the field at which the pivot mode leaves the horizon. 
It can be rescaled by a shift 0» — > 0» + and the range over 
which changes during inflation varies greatly between mod- 
els. Consequently we treat the remaining number of e-folds, A», 
after the pivot scale leaves the horizon as a free parameter with 
a wide uniform prior, since this quantity has a consistent inter- 
pretation across models. The pivot scale used to compute A, is 
K - 0.05 Mpc However, given our ignorance concerning the 
epoch of entropy generation, a multitude of entropy generation 



scenarios for each inflationary potential can occur Some possi- 
bilities are as follows, with parameters referring to Eq. 24. 

1. Instantaneous entropy generation. 

At the end of inflation, all the energy in the inflaton field is 
instantaneously converted into radiation. 

2. Restrictive entropy generation scenario (narrow range for 

Wint)- 

p'/"* = 10'' GeV, and Wi„t e [-1/3, 1/3]. 

3. Permissive entropy generation scenario (wide range for 

Wint). 

= 10^ GeV, and Wi„t e [-1/3, 1]. 

The equations of state with Wint in the range [1/3, 1] appear less 
plausible, but models with these values have been put forward 
(Pallis, 2006), so this cannot be completely excluded. Moreover 
the Wint parameterization captures a variety of scenarios in which 
the post-inflationary Universe is thermalized, but not radiation- 
dominated, including phases of coherent oscillations (Martin & 
Ringeval, 2010; Easther et al., 2011), resonance (Traschen & 
Brandenberger, 1990; Kofman et al., 1997; Allahverdi et al., 
2010), kination (Spokoiny, 1993; Chung et al., 2007), secondary 
or thermal inflation (Lyth & Stewart, 1996), moduli domination 
(Banks et al., 1994; de Carlos et al., 1993), primordial black hole 
domination (Anantua et al., 2009), or a frustrated cosmic string 
network (Burgess et al., 2005), all of which lead to an expansion 
rate differing from that of a radiation-dominated universe. 

At the other extreme, the decision to exclude Wint < -1/3, 
as done here, is not completely justifiable. We cannot, for ex- 
ample, rule out a first-order phase transition at a lower energy 
scale that would drive vvint below -1/3, but here we neglect this 
possibility. Our analysis does not preclude a secondary period 
of inflation, but does require that the average expansion during 
the post-inflationary regime parameterized by Wim should not be 
inflationary. This caveat should be kept in mind. 

For some of the parameterized models, tighter constraints 
can, in principle, be placed on wim. It has been argued (see e.g., 
Liddle & Leach, 2003) that for the A(p'^ potential, the uncertain- 
ties concerning entropy generation contribute almost no uncer- 
tainty in the determination of 0*. This is because according to 
the virial theorem, a field sloshing about = in a quartic po- 
tential has the same average w, namely w - 1/3, as the radiation 
equation of state. More generally, for a potential of the form 0" 
around the minimum, vvvii- = (w) = (n - 2)/{n + 2) (Turner, 
1983); therefore, one may argue that wim should be restricted 
to the interval whose endpoints are 1 /3 and Wvii . This approach 
was taken by Martin & Ringeval (2010) in obtaining Bayesian 
constraints on the reheating temperature for monomial poten- 
tials from the CMB. However this scenario requires a carefully 
tuned potential which has approximately a 0" shape both at large 
field values, and near the origin, far below the inflationary scale. 
Typically, potentials for which V((f)) ~ (p" at large field values 
can have very different shapes near the origin. Thus, following 
Easther & Peiris (2012), in this paper we choose to explore a 
broader range of Wint for these models (including the cases above 
as subsets), in order to obtain data-driven constraints on Wi„t. 

In this paper we focus on the three representative scenarios 
itemized above, referred to hereafter as: (1), (2), and (3). Our 
algorithm draws a value of A* and then, given the value of pth, 
computes Wim, which is a derived parameter — models for which 
Wint lies outside the specified range of each scenario under con- 
sideration are excluded. 

Our full set of priors for the inflationary physics is given in 
Table 6. Dimensionful quantities are expressed in units with re- 
duced Planck mass Mpi set to unity. 
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Model Instantaneous Restrictive Permissive 

entropy generation entropy generation entropy generation 





ln[fi/fio] 




ln[£/£o] 




ln[£/fio] 




n = 4 


-14.9 


25.9 


-18.8 


27.2 


-13.2 


17.4 


n = 2 


-4.7 


5.4 


-7.3 


6.3 


-6.2 


5.0 


n = 1 


-4.1 


3.3 


-5.4 


2.8 


-4.9 


2.1 


n = 2/3 


-4.7 


5.1 


-5.2 


3.1 


-5.2 


2.3 


Natural 


-6.6 


5.2 


-8.9 


5.5 


-8.2 


5.0 


Hilltop 


-7.1 


6.1 


-9.1 


7.1 


-6.6 


2.4 


ACDM 


-4940.7 


9808.4 











Table 7. Inflationary model comparison results. For each model and set of assumptions concerning entropy generation [(1), (2), (3)], 
the natural logarithm of the Bayesian evidence ratio as well as ^xlf{ for the best-fit model in each category are indicated, relative to 
the ACDM concordance model (denoted by subscript "0"); In &o -2 In Xo for the latter are also given. 



Due to the non-trivial likelihood surfaces and the large di- 
mensionality of the parameter spaces explored in this section, 
we use ModeCode coupled with MultiNest v3 . 0^ to map 
out the parameter space. In addition to the standard nested sam- 
pling (NS) algorithm. Mult iNest v3 . enables nested impor- 
tance sampling (NIS), resulting in substantial speed gains** and 
significantly enhancing the accuracy of the Bayesian evidence 
computation over NS alone for the same computational setup. 

5. 1. Results 

Table 7 presents model comparison results for the ensemble 
of parameterized potential families described above. We report 
the Bayesian evidence (model-averaged likelihood) ratio, which 
provides a self-consistent framework for calculating the betting 
odds between models (see Sect. 3.4). The uncertainty in these 
log evidence values is ^ 0.2. We also report the A^f^^. values 
computed from the 2 In -Cmax values found by the sampler 

The monomial models have a single-parameter potential, and 
the natural and hilltop inflation models have two parameters 
each. All reheating scenarios except case (1) contribute one ad- 
ditional parameter to the inflationary sector. The evidence ratios 
and ^xlf[ values are presented with respect to the ACDM cos- 
mological model. 

None of the inflationary models tested here fit the data as 
well as the ACDM model. This mostly reflects that there is 
no evidence in the data for r different from zero. Furthermore, 
the priors listed in Table 6 for the ACDM primordial sector 
are purely phenomenological, roughly corresponding to ranges 
somewhat broader than WMAP constraints. Nan^owing them 
around the best-fit model arbitrarily increases the evidence in its 
favour. Instead it is instructive to compare the relative evidence 
for the inflationary models presented. 

Table 7 shows that the Acp'^ model is decisively ruled out by 
Planck, confirming previous analyses by the WMAP team (Peiris 
et al., 2003; Spergel et al., 2007; Dunkley et al., 2009; Komatsu 
et al., 201 1) based on the model track plotted on the n^-r plane. 
Recent model selection analyses (Martin et al., 2011; Easther & 
Peiris, 2012) with WMAP7 found that the model was akeady 
under severe pressure, disfavoured by odds of ~ 400: 1 against. 



' Made available ahead of public release to the Planck Collaboration 
by Farhan Feroz and Mike Hobson. 

We have made extensive tests of NIS vs. NS, and chosen the fol- 
lowing settings for the computations presented here: NIS on, constant 
efficiency mode on, 300 live points, tolerance and efficiency parameters 
set to 0.5 and 0.02, respectively. 



With Planck, the odds against this model are at least 500, 000: 1 
when compared with ACDM, for a very broad set of reheating 
scenarios. The same conclusion is reflected by the extremely 
poor ^xlff values for the model. Given the strength of our re- 
sults in the flexible setting of the permissive entropy generation 
scenario, it is possible not just to rule out models where the po- 
tential is of the quartic form in the full range from the origin to 
the inflationary scales, but also a general class where the poten- 
tial is of the n - A form in the 0-range where the cosmological 
perturbations are generated, but exhibits a different shape near 
the origin. 

Two other large field models, the quadratic potential and nat- 
ural inflation, are under some pressure from Planck data, espe- 
cially when broader entropy generation scenarios are considered. 
Compared with the ACDM model, these models are disfavoured 
by ^^x\s ~ 5-6 depending on the reheating scenario. This re- 
flects the analysis of Sect. 4, where the overlap of the model 
predictions and the data constraints on the n^- r plane is seen to 
be mostly outside the joint 68% CL contour. However, from the 
Bayesian evidence point of view, it is too early to declare these 
models incompatible with the data. To make this judgement, it 
is more conservative to compare these models with the n - \ 
case, which has the best evidence with respect to ACDM, rather 
than with ACDM itself, which provides our reference point for 
the evidence calculation, but has arbitrary prior ranges. In their 
simplest forms - instant entropy generation - the n - 2 and nat- 
ural inflation models are only disfavoured by odds of ~ 1-12:1 
against, which does not rise to a high level of significance''. 

The models most compatible with the Planck data in the set 
considered here are the two interesting axion monodromy poten- 
tials, n - 1 (McAllister et al., 2010), and n - 2/3 (Silverstein 
& Westphal, 2008), which are motivated by inflationary model 
building in the context of string theory. The /? = 4 hilltop model 
presents an interesting case. This model was previously found to 
be compatible with WMAPl data, performing almost as well as 
the monodromy potentials (Easther & Peiris, 2012). However it 
exhibits significant tension with the Planck data, both in terms 
of evidence ratios and the maximum likelihood. The only ex- 
ception is the entropy generation scenario (3) which has odds of 
> 1000:1 against compared to ACDM, and yet the maximum 
likelihood is not significantly different to the n - \ case. This 
indicates that while the extra freedom allowed by the least con- 
servative reheating scenario improves the best fit, this prior is 
not very predictive of the data. However, the result is counterin- 



' In comparison, odds of 150:1 are considered highly significant in 
this context. 
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Fig. 8. Constraints on «s vs. r at kt - 0.002 Mpc ' for the inflationary models considered (i.e., power law potentials with 
n - 2/3, 1,2, and 4, natural inflation, and hilltop inflation), showing joint 68% and 95% CL. Blue and grey distributions corre- 
spond to the restrictive and permissive entropy generation scenarios, respectively. The instant reheating case corresponds to the thin 
(red) contours in the natural and hilltop panels; for the single parameter models, this case corresponds to the lowest-r extremity of 
restrictive case. The difference between the natural inflation region on Fig. 1 and the natural inflation constraints shown here is due 
to the strong projection effect described in the text. 




Fig. 9. Constraints on «s vs. Wmt at - 0.002 Mpc ' for the inflationary models considered, as in Fig. 8. The instant reheating case 
(1) corresponds to Wint = 1/3. 

tuitive and merits further comment; we will consider it further We now turn to parameter constraints. Fig. 8 presents 
at the end of this section when parameter estimation results are marginalized joint constraints from Planck alone on the derived 
discussed. parameters ns,o.oo2 and ro.002; Fig. 9 shows the corresponding 

joint constraints on ^5,0.002 and Wint, again a derived parameter. 
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Fig. 10. Constraints on ln(A) vs. Wint for the inflationary models considered, as in Figs. 8 and 9. 



and Fig. 10 shows joint constraints on Wint vs. the potential pa- 
rameter log(/l) (log(A) in the case of natural inflation). It is in- 
structive to consider the three sets of figures together. The re- 
strictive and permissive entropy generation scenarios are shown 
on all panels; the instant reheating case is shown for the two pa- 
rameter models, natural inflation and hilltop, in Fig. 8 - for the 
monomial potentials this case corresponds to the lowest-r ex- 
tremity of the restrictive case constraints, and to vvint - 1 /3 in 
the other two figures. 

The quartic potential conflicts with the data because it pre- 
dicts a high tensor-to-scalar ratio. Hence the model maximizes 
its likelihood by pushing towards the lower-r, bluer-Wj limits 
of its parameter space, which corresponds to increasing Wint as 
allowed by the reheating prior. The contours terminate at the 
lowest-r limit when each reheating case hits its vvint prior upper 
limit (i.e., vvint = 1/3 and Wint = 1 for restrictive and permis- 
sive entropy generation, respectively). For each case, the lower 
limit on Wint, corresponding to the reddest-Ws extremity of the 
confidence contours, is data-driven. The quadratic potential en- 
counters the same difficulty but at a less extreme level. 

The two axion monodromy potentials are compatible with a 
wide range of entropy generation scenarios. The instantaneous 
reheating scenario is compatible with the data for both models. 
For restrictive entropy generation in the « = 1 model, we obtain a 
data-driven upper limit on Wint, which just touches the Wint = 1/3 
case at the 95% CL. At the lower limit, the Wint posterior is trun- 
cated by the prior, as for the n - 213 case. For the latter, there 
is a data-driven upper limit on Wiai which is controlled by the 
upper limit on n^. For permissive entropy generation, the upper 
and lower limits on w^^i for both models are data-driven, corre- 
sponding to the upper and lower limits on respectively. 

The constraints on natural inflation require some interpreta- 
tion. The relationship between the empirical and r parameters 
and the potential parameters for natural inflation is discussed in 
detail by Savage et al. (2006) and Mortonson et al. (201 1) along 
with parameter constraints derived from WMAP3 (Spergel et al., 
2007) and WMAPl (Larson et al., 2011), respectively. In this 



model, there is a degeneracy between / and A in the limit where 
these parameters are large so that natural inflation resembles the 
quadratic model. Our priors are chosen to exclude most of this 
region. The priors on log(/) and log(A) still allow a region of 
nearly degenerate models that contribute to a ridge seen in the 
natural inflation panel of Fig. 8. These models closely match the 
values of and r seen in the quadratic potential constraints. The 
marginalized constraints on and r depend strongly on the prior 
on log(/) due to the projection of a large number of degenerate 
models into this ridge. Therefore the apparent preference for this 
region of parameter space over models with lower values of r 
is largely due to this effect and is not driven by the data. This 
highly nonlinear mapping between the logarithmic priors on the 
potential parameters and the power law parameters (which are 
derived parameters in this analysis) leads to a strong projection 
effect that accounts for the difference in visual appearance be- 
tween these contours and the region labeled 'natural inflation' in 
Fig. 1. 

Generally, for fixed Nt, decreasing A and / reduces both 
and r. Thus natural inflation models can have lower values of r 
than the quadratic potential without increasing Ws and A^,. This 
feature means that the potential parameters for this model are 
relatively uncorrected with Wint, in contrast with the other mod- 
els considered here, as illustrated in the Fig. 10. Nevertheless 
we obtain data-driven bounds on in the permissive reheating 
case as well as a lower bound in the restrictive reheating case. 
Both bounds overlap with the instantaneous reheating limit. 

The left panel of Fig. 1 1 shows our lower limit on / > 10.0 
Mpi (95% CL), compared with the WMAPl Umit / > 5.0 Mpi 
(95% CL) reported by Mortonson et al. (2011). Indeed, the 
Planck limit is in agreement with the Planck prediction pre- 
sented in that work. There is a hint of an upper limit on / as 
well, driven by the fact that this corresponds to the quadratic in- 
flation limit, which is in tension with the data. However this is 
only a 1 cr effect. 

The p - A hilltop model has two distinct branches: a small- 
field scenario, where r < 0.001 and < 0.95, and a large- 
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Fig. 11. Potential parameters for natural inflation and hilltop in- 
flation, as in Figs. 8, 9, and 10. On the natural inflation panel, 
instant reheating corresponds to a thin diagonal along the top 
edge of reheating case (3). 



^eW limit in which V(<p) ~ 4> (Adshead et al., 2011). Physically, 
the small-field limit is consistent with the Lyth bound, and we 
can select it by fixing log(A) < -2.5 in the prior We observe 
that the data select out this small field branch, which requires 
explanation given that we know that the n - 1 model (the limit 
of the large field branch) is perfectly compatible with the data. 
In fact, this can be tested by restricting the log(A) prior by hand 
to the large-field branch. The hilltop model constitutes a difficult 
sampling problem, as is apparent from Fig. 4 in Easther & Peiris 
(2012). An examination of the nested sampler's progress in the 
case of the full hilltop prior reveals the reason for the small- 
field branch being selected out. In this case, the posterior for the 
large-field branch is extremely thin compared to the model prior 
in this regime - much thinner than the posterior for the small- 
field branch in comparison to its respective prior. Therefore this 
region occupies very little probability mass, and is dropped in 
preference to the more predictive small-field branch. This high- 
likelihood but extremely thin ridge is also responsible for the 
counterintuitive result reported in the model-selection analysis, 
where the model was found to have a good ^xl^f with respect to 
ACDM, and yet be highly disfavoured by the Bayesian evidence. 



In summary, confirming the results of Sect. 4, this model is in 
agreement with the data in the limit where it overlaps the linear 
« = 1 model. But since this region of high likelihood occupies a 
very small fraction of the prior, this model is heavily penalized 
by the Bayesian evidence for failing to predict the data over most 
of its prior space. 



6. Observable window of inflation 

Section 4.2 presented an analysis of several representative in- 
flationary potentials within the framework of the slow-roll ap- 
proximation, and their compatibility with the Planck data. The 
results are summarized in Fig. 1. In that case the full potential is 
considered in order to identify a plausible range for the location 
of 0t on the potential V{(p). This requires a complete story: i.e., 
the potential must be specified starting above the point where 
the largest observable scales first exited the Hubble horizon, and 
extending to the minimum of the potential. 

In this section we explore another approach. We adopt the 
point of view that we are interested in reconstructing the infla- 
tionary potential only over the observable range — that is, the in- 
terval of <p corresponding to the scales observable today in the 
CMB. We constrain the potential over the range where these 
scales exited the horizon during inflation, as well as a few e- 
folds before and after - the cosmological perturbations are not 
imprinted instantaneously at the moment of horizon crossing, 
but rather, gradually over a few e-folds. We expand around (pt, 
taking the view that a plausible extension of the potential out- 
side this observable range is always possible, so that one has 
precisely the number of e-folds of inflation needed for 0» to cor- 
respond to kt today. 

The argument is that one can always end inflation abruptly 
by imposing a sharp waterfall feature where needed, or prolong 
inflation by inserting a sufficiently long plateau into the poten- 
tial by hand, e.g., for models with a large tilt. A foreseeable ob- 
jection to this approach is that the extensions of the potential 
required outside the observable window may render the poten- 
tial unnatural. This possibility should be kept in mind, although 
naturalness is an elusive and uncomfortably subjective concept. 
The analysis in this section does not rely on the slow -roll approx- 
imation. Instead each k mode is integrated exactly by numerical 
integration, as described in Sect. 2.2, under the assumption of a 
canonical kinetic term. 

Two complementary approaches have been taken in the lit- 
erature to reconstruct the potential. In the first approach, which 
we follow in this paper, we expand the potential V{^) directly in 
powers of (0 - 0*). In this case the numerical integration of the 
slow-roll solution must start sufficiently early so that any initial 
transient has had a chance to decay, and one is in the attractor so- 
lution when the dynamics of the largest observable modes in the 
Universe today start to have an interesting evolution. This is the 
approach followed, e.g., by Lesgourgues & Valkenburg (2007); 
Mortonson et al. (201 1) using publicly available codes'". 

A second approach expands H{(fi) as a Taylor series in {(p - 
(pt). As discussed in Sect. 2.2, this has the advantage that i/(0) 
determines both the potential y(0) and the solution (pit), so the 
issue of having to start sufficiently early in order to allow the ini- 
tial transient to decay is avoided. This method was used in, e.g., 
Kinney (2002), Kinney et al. (2006), Peiris & Easther (2006a), 
Easther & Peiris (2006), Peiris & Easther (2006b) and Peiris 



http : / /wwwlapp . in2p3 . f r /valkenbu/inf lationH/ , 
http : //zu server 2 .star.ucl.ac. uk/ ~ hi r any a /Mode Code/ . 
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Fig. 12. Posterior distribution for the first three slow-roll poten- 
tial parameters using Planck+WP data. In the n = 2-4 cases, the 
inflaton potential is expanded to the «th order, and the spectrum 
is obtained by fitting the numerically-computed power spectrum 
to the data, with no slow-roll approximation, and no assumption 
about the extrapolation of the potential outside of the observable 
window. 
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Fig. 14. Observable range of the best-fitting inflaton potentials, 
when V((f>) is Taylor expanded at the «th order around the pivot 
value (pt, in natural units (where y/SnM^i - 1), assuming a 
flat prior on ey, rjv, ^y, and Wy, and using Planck+WP data. 
Potentials obtained under the transformation ((p-<pt) — > (0, - <p) 
leave the same observable signature and are also allowed. The 
sparsity of potentials with a small Vo - V((f>f) comes from the 
flat prior on £y rather than on ln(y()): in fact Vo is unbounded 
from below. 

Table 8. Numerical reconstruction of potential parameters, com- 
pared to results with the slow-roll approximation, when ten- 
sors and running are included (Planck+WP 95% CL, with k^ = 
0.05 Mpc '). The effective value is given relative to the model 
with a quadratic potential. 
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Fig. 13. Posterior distribution for the Taylor expansion coeffi- 
cients, y,, of the inflaton potential. The potential is expanded at 
nth order, assuming a flat prior on ey, rjy, ^y, ml. The coeffi- 
cients y, are expressed in natural units (where y/SnM^i = 1). 
The contours show only half of the allowed regions for potential 
parameters: the other half is symmetric, with opposite signs for 
Vi and V3. 



& Easther (2008), using analytic and semi-analytic approxima- 
tions, and in Lesgourgues et al. (2008), Powell & Kinney (2007), 
Hamann et al. (2008c) and Norena et al. (2012), using a fully nu- 
merical approach. 

These approaches could lead to results that differ from those 
in Sect. 4.2. Firstly, if the running of the index is large, the slow- 
roll approximation taken to second order is not necessarily accu- 
rate for all models allowed by the data. The relation between the 
spectral parameters (In As, «s, dn^/dlnk, r) and the underlying 
inflationary potential y(0) is therefore uncertain. Secondly, for 



spectra with a large running, there is no guarantee that an infla- 
tionary model giving such a spectrum exists. All allowed mod- 
els have r <K 1, so these models are consistent with ey{kt) <K 1 
at the pivot scale. However, towards the edge of the observable 
range, the potential may become incompatible with €y{k) < 1 
(i.e., with the requirement of inflationary expansion). These pos- 
sible pitfalls are avoided using the methods in this section, since 
the data have been fit directly by the candidate y(0) or H((p), 
computed numerically without any slow roll assumptions over 
the entire observable range. 

We define a class of models over the observable range based 
on the expectation that the potential should be smooth. y(0) is 
approximated by a Taylor expansion up to order n, and we ex- 
plore the cases n = 2, 3, and 4. For each V{4>), we integrate over 
inflationary fluctuations using the inflation module implemented 
in class" (Lesgourgues, 201 1; Bias et al., 201 1), as described 
in Lesgourgues & VaUcenburg (2007). Potentials are rejected for 
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Table 9. The scalar amplitude, tilt, running, running of the run- 
ning, and tensor-to-scalar ratio inferred from a numerical recon- 
struction of the inflaton potential (Planck+WP 95% CL, with 
= 0.05 Mpc-i). 
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which the attractor solution cannot be reached when the largest 
observable scales cross the Hubble radius. The parameters we 
sample are the potential and its derivatives at the pivot scale, 
when kt crosses the Hubble radius during inflation. 

To avoid parameter degeneracies, we impose uniform pri- 
ors on ey, riv^^y, and -uTy at the pivot field value The advan- 
tage of uniform priors on these parameters is that - to the extent 
that the slow-roll conditions are satisfied - these coefficients re- 
late linearly to observable quantities such as «s, r, dnj/dlnA: and 
d^Hs/dlnA;^. Fig. 12 and Table 8 show the posterior probability 
for these coefficients, and Fig. 13 shows the posterior probabil- 
ity for the Taylor series coefficients V,-. In Fig. 14, we show the 
observable range of the best-fitting inflaton potentials (for a sam- 
ple extracted randomly from the converged Markov chains). The 
edges of the observable range correspond to Hubble crossing for 
the minimum and maximum values of k used in the Boltzmann 
code. 

When fitting V{(p) for each model in parameter space, we 
compute (hs, dns,/d\nk, cfiiis/dlnk^, r) at the pivot scale a pos- 
teriori directly from the numerical primordial spectra. The re- 
sults are shown in Table 9 and can be compared with those of 
Table 5 in Sect. 4. We can also use the results of Sect. 4 for the 
ACDM-HrH-dns /d In model to infer the potential parameters (ey, 
Tjv, ^y) using the second-order slow-roll expressions, and com- 
pare different approaches in the space of potential parameters 
(see Fig. 12 and the last column in Table 8). 

The model with a quadratic potential in the observable win- 
dow (n - 2) leads to bounds on ey, rjy, «s, and r very close to 
the ACDM-Hr case. This is not a surprise since such potentials 
cannot simultaneously give values of «s and r compatible with 
the data, and a large running. A significant dws/d In k can be gen- 
erated only in the presence of a large (i.e., with a significant 
V"). Since quadratic potentials produce little running, they are 
faithfully described by the slow-roll approximation. 

The model with a cubic term (n = 3) has the freedom 
to generate a large running, dn^/dlnk. Indeed one can check 
that the results for the « = 3 model are close to those of the 
ACDM-HrH-dns/d In model presented in Sect. 4. The agreement 
between these two models remains very good, despite the fact 
that in the presence of a large running, the slow-roll approxi- 
mation can become inaccurate. The running in a potential with 
n = 3 is not exactly scale invariant; this is not captured by the 
ACDM-HrH-dns /dink parameterization. 

The n - 4 model has even more freedom, allowing a 
considerable running of the running d^n^/dlnk^ (to the extent 
that inflation holds during the observable e-folds). In that case, 
the spectrum is better fitted when the two parameters r and 
d^«s/dlnA:^ are non-zero. In Fig. 14, we see that most n - A 



potentials have a long and steep tail for (p < cf>t, with a kink 
around 0* - 0.4 (in natural units). This shape generates a signif- 
icant running on the largest observable scales, while preserving 
a smaller running on smaller scales. With such a feature in the 
scalar primordial spectrum at large scales, combined with a non- 
zero contribution from tensor fluctuations, the best-fit model for 
n - 4 has a temperature spectrum very close to that of the min- 
imal ACDM model for £ > 40, but not for smaller multipoles: 
the amplitude of the Sachs-Wolfe plateau is smaller. This allows 
the large-scale data points from Planck to be fit slightly better. 
However, the case r - dns/dlnA: - still lies at the edge of the 
95% CL, and the minimum effective of this model is smaller 
than in the n - 2 case by only 3.7. 

A comparison of the n - 3 and n - A results clearly shows 
that the process of expanding the inflation potential at various 
orders and fitting it to the data does not converge (at least not by 
n - 4). Given the 1 - 2cr preference of Planck data for a non- 
zero running and running of the running, we find that a model- 
independent reconstruction of the inflaton potential is not pos- 
sible under the assumptions of this section. In other words, as 
long as we assume that V{<p) can be described during and after 
observable inflation by a polynomial of order 2 or 3, we can put 
strong bounds on ey, rjy. But if we introduce more derivatives to 
describe the observable part of the potential, and keep complete 
freedom to extrapolate V{(p) outside this region, the constraints 
can be easily evaded. 

7. Primordial power spectrum reconstruction 

In this section we report on a search for features in the pri- 
mordial power spectrum. In the basic six parameter model 
studied in the companion Planck paper Planck Collaboration 
XVI (2013), the primordial power spectrum PKik), which 
includes only the adiabatic mode, is modeled using the 
power law Vnik) - A^iklkt)"'^"^ , for which the best fit val- 
ues are As - 2.20 x 10"' and = 0.9603 for a pivot scale 
A;» = 0.05 Mpc"'. An extension of this parameterization is 
also considered allowing for a running of the spectral index 
(dns/dlnA; = -0.013 + 0.009). But in all cases considered it was 
assumed that the power spectrum is smooth and without bumps, 
sharp features, or wiggles. In this section we investigate whether 
any statistically significant evidence for features is present in the 
data when these assumptions are relaxed. Allowing an arbitrary 
function for the input power spectrum is not an option because 
in this case the recovered primordial power spectrum is domi- 
nated by small scale noise. Instead we consider here a penalized 
likelihood approach where a preference for smooth power spec- 
tra is imposed. The companion Sect. 8 of this paper pursues a 
complementary approach where several parametric models for 
wiggles and features are explored to see whether a statistically 
significantly better fit can be obtained. 

An extensive literature exists on how to search for features 
in the power spectrum using a wide range of methods. The fol- 
lowing papers and the references therein provide a sampling 
of the literature on non-parametric reconstruction: Richardson- 
Lucy deconvolution (Lucy, 1974; Richardson, 1972; Hamann 
et al., 2010; Shafieloo & Souradeep, 2004, 2008), deconvolution 
(Tocchini-Valentini et al., 2005, 2006; Ichiki & Nagata, 2009; 
Nagata & Yokoyama, 2008, 2009), smoothing splines (Verde & 
Peiris, 2008; Peiris & Verde, 2010; Sealfon et al., 2005; Gauthier 
& Bucher, 2012), linear interpolation (Hannestad, 2004; Bridle 
et al., 2003), and Bayesian model selection (Bridges et al., 
2009; Vazquez et al., 2012). The approach pursued here follows 
Tocchini-Valentini et al. (2006) and Gauthier & Bucher (2012) 
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Fig. 15. Planck primordial power spectrum feature search results. Top four panels: The reconstructed power spectrum at four values 
for the smoothing parameter A. The red curves indicate the maximum likelihood configuration for the fractional deviation f(k) of 
the power spectrum relative to a power law fiducial model (with = 2.20 x 10"^ and = 0.9603) for the penalized likelihood. 
The error bars have a width corresponding to the minimum reconstructible width (the minimum width for a Gaussian feature so that 
the mean square deviation of the expectation value of the reconstruction differs by less than 10%), and a vertical extent showing the 
Icr and 2cr limits for the fractional deviation averaged over the box. The grey hashed regions at the far left and right show where 
the fixing prior (i.e., a) sets f{k) - 0. The inner grey regions show where the reconstruction bias is so great that the minimum 
reconstructible width is undefined. Mock features in this region have reverberations over the entire interval. With A - and 
A - 10^, we find statistically significant fluctuations around k ~ Q. \ Mpc"'. Lower panels: The lo" error bars for three combinations 
of cosmological parameters at the four values of A. The maximum likelihood value for the fiducial model is indicated by the dashed 
line for comparison. 
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most closely. More technical details and extensive tests validat- 
ing the method can be found in the latter reference. 

Let Poik) = As(k/kt)"^^^ be the best fit power spectrum of the 
six parameter model. We define a general ansatz for the power 
spectrum in terms of a fractional variation, f{k), relative to this 
fiducial model, so that 

r^^mkyli+m]. (4?) 

Any features are then described in terms of f{k). 

In this analysis we use the Planck+'WP likelihood supple- 
mented by the following prior, which is added to -2 In X.- 

dKf(K) + a I dKf(K). 

where k = Ink. The first term penalizes any deviation from a 
straight line of the function /(k), and the second term fixes the 
/(k) to zero where there are effectively no constraints from the 
data. The value for A controls the smoothness of the reconstruc- 
tion, but the precise value of a is less important. It must be large 
enough to force f{k) towards zero when k < K,^in and k > /fmax 
but no so large as to render the matrices ill-conditioned. We use 
a = 10\ 

We represent f(k) using a cubic B-spline on a grid of points 
in ^-space uniformly spaced in a- = In A: with step size A/c - 0.025 
and extended from k = -12.5 to a: = -0.3 giving us a total of 485 
knots so that f = {./i lfJi • The density of grid points is sufficiently 
large so that artefacts around the scale of the knot spacing are 
suppressed for the values of A used here. Given the large number 
of dimensions, it is not practical to explore the likelihood using 
MCMC methods. However, for the power spectrum parameters, 
the predicted CfS are related by a Unear transformation given 
fixed cosmological parameters, allowing us (for fixed cosmolog- 
ical and nuisance parameters) to find the maximum likelihood 
solution using the Newton-Raphson method.'^ We define 

M(0)= min {-2ln j:(Q,t)+t'^R(A,a)t\ (49) 

where the vector represents the non-power spectrum cosmo- 
logical parameters and foreground nuisance parameters, and we 
first minimize over / using Newton-Raphson iteration and then 
in an outer loop minimize over using the downhill simplex 
algorithm. To achieve this, the Planck likelihood code was mod- 
ified to compute the gradient and Hessian of the likelihood with 
respect to the CfS. 

The cosmological Boltzmann solver CAMB was modified to 
accept the vector of primordial power spectrum knots f. By de- 
fault CAMB calculates the CfS for a subset of t and interpolates 
to obtain the full multipole power spectrum. Instead we calculate 
CfS at each { explicitly. 

The boundaries K,^i„ and /Cmax defining where /{k) is allowed 
to differ from zero are chosen to match the range of { con- 
strained by the high-/" likelihood. The likelihood includes C^s 
between i = 50 to ^ = 2500, which roughly corresponds to 
k € [0.003,0.2] Mpc^'. The low-{ likelihood covers ^ = 2 to 
£ = 49, which roughly corresponds to k e [10 0.003] Mpc 

Since the Planck likelihood uses a quadratic approximation, the 
maximum likelihood solution can be found by solving a linear system. 
Newton-Raphson, however, is needed for other CMB likelihoods which 
include terms beyond quadratic order. 
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Fig. 17. CMB multipole spectrum residual for the primordial 
power spectrum test feature. The test feature (top) is set to the 
anomalously large deviation of the primordial power spectrum 
reconstruction for A - 10^ in the interval 0.1 Mpc ' < k < 
0.15 Mpc and is zero elsewhere. Bottom: The angular spec- 
trum corresponding to this feature. We observe a large dip at 
£ ^ 1800. 

In this range of {, cosmic variance is large, making feature 
detection difficult. Calculating the gradient and Hessian of a 
pixel-based low-^ likelihood is computationally time consum- 
ing. We therefore use the low-^ likelihood only to constrain the 
cosmological parameters. We choose kj„i,^ - 0.005 Mpc"' and 
^max = 0.3 Mpc"'. Within this k range, variations in the C(S due 
to the fi are too small to affect the overall likelihood through the 
low-^" likelihood. We observed that the difference in the low-^ 
likelihood between the reconstructed f{k) and f{k) = is small 
(< 1%) compared to the difference in the high-^ likelihood. 

The cosmological parameters t and are almost completely 
degenerate for the temperature anisotropy, except at very low {, 
so we fix T to its best fit value for the fiducial model. The likeli- 
hood contains additional nuisance parameters that model fore- 
ground components and beam shapes, as discussed in Planck 
Collaboration XV (2013). Many of the nuisance parameters, it 
can be argued, are unlikely to introduce spurious small scale 
structure because they represent foreground models with a power 
law and thus smooth angular power spectrum. However some 
nuisance parameters, in particular those describing beam un- 
certainties, could conceivably introduce artefacts into the re- 
construction. Unfortunately, converging to the correct maximum 
likelihood reconstruction with all the beam shape parameters in- 
cluded is prohibitively time consuming. Therefore we fix the 
nuisance parameters to their fiducial best fit values, leaving a 
more detailed examination of this issue to future work. 

We found that simultaneously allowing extra degrees of free- 
dom for small scale structure and variations in the cosmological 
parameters changes the best fit fiducial model — that is, and 
«s — so that the variations with respect to the fiducial model no 
longer visibly gave the best straight line fit. Therefore we allow 
As and to vary, so that the fiducial model is indeed the best 
straight line fit through the plotted data points. Detailed investi- 
gation showed that neither the priors, nor low-^, nor high-/' data 
play an significant role in determining the best fit fiducial model. 
This effect is small and within the error bars for As and estab- 
lished assuming the fiducial model. To summarize, we maximize 
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Fig. 16. CMB multipole spectrum residuals for best fit primordial power spectrum reconstruction at A - 10^. The panels show the 
Cf spectrum residuals (compared to the best fit power law fiducial model represented by the horizontal straight dashed line) for the 
four auto- and cross-spectra included in the high-/" likelihood. Here, Dg = {(({ + l)/(2:7r))Cf. The data points have been binned with 
A{ = 31 and foregrounds subtracted according to the best fit foreground parameters. The solid black line shows the CMB spectrum 
residual for the maximum likelihood primordial power spectrum reconstruction with smoothing parameter A - 10^. 



the likelihood with respect to the control points // and the three 
cosmological parameters h, Q.^h-, and Q.\jh^. We update the fidu- 
cial model (As and «s) at each iteration by finding the best fit 
power law through the current best fit reconstruction. 

Once the maximum likelihood solution has been found, the 
second derivatives about this solution are readily calculated by 
extracting the relevant matrices for most of the components, and 
estimating the remaining components using finite differences. 
The second derivative matrix is used to estimate the error on 
the reconstructed fi and the three cosmological parameters h, 
Q.^h^, and Qb/z^. Monte Carlo simulations of a fiducial data set 
with a simplified CMB likelihood including some of the non- 
Gaussianities ' suggest that this approximation of the error is 
accurate for A>\(fi . 

Fig. 15 summarizes our results, showing the estimated f{k) 
in bins, and corresponding Icr and 2cr errors. Errors in k are also 
shown to represent the minimum reconstructible width evalu- 
ated at the middle of each bin. This is the minimum width that 
a Gaussian feature must have to be reconstructed with a small 
enough bias such that the mean reconstruction differs by less 
than 10% rms. The minimum reconstructible width is closely 



" The simplified likelihood -llnLt^^ = Yl^^C^t + 1)^^^ - 
In — 1 ) assumes full-sky coverage and isotropic instrument noise. 



related to the correlation length, so that the errors between ad- 
jacent bins are weakly correlated and the total number of bins 
represents roughly the effective number of independent degrees 
of freedom. 

While the plots with a significant roughness penalty — that 
is, with A - 10^ and A - 10^ — do not show any statistically 
significant evidence for features standing out from the noise of 
the reconstruction, for a smaller roughness penalty — that is, for 
A - 10"* and A - 10^ — a nominally statistically significant fea- 
ture is clearly visible around k x 0.13 Mpc"'. We do not un- 
derstand the origin of this feature, which may be primordial or 
may arise as a foreground or other systematic error in the high-/' 
portion of the likelihood. It should be noted that most of the ro- 
bustness tests described in the likelihood paper assume smooth 
power spectra. The maximum excursions are locally at 3.2cr and 
3.9cr for A - 10^ and A = 10^, respectively. In each of these two 
cases we correct for the look elsewhere effect by calculating the 
probability that one of the plotted error bars deviates by the same 
number of or more standard deviations. This calculation is car- 
ried out using the covariance matrix of the plotted error bars. We 
obtain p = 1.74% and p - 0.21%, which corresponds to 2.4cr 
and 3. Icr, respectively. Additional simulations were carried out 
to validate the method by generating mock data according to the 
fiducial model and measuring the errors of the reconstruction 
obtained. These investigations confirm the error model. These 
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tests where carried out both with and without test features. It can 
be argued that foregrounds are unUkely to explain the observed 
feature because all the foreground models involve smooth power 
law templates whereas this feature is localized in multipole num- 
ber. It is important to assess by means of a more extensive set of 
simulations whether the statistical significance assigned to this 
result is accurate. 

We investigate which CMB angular multipoles correspond 
to this apparent feature. Fig. 16 shows the Q residual from 
the reconstructed power spectrum with the best fit power law 
power spectrum subtracted together with the data for each fre- 
quency map correlation combinations used in the CamSpec 
likelihood. We observe a smooth dip around I ^ 1800, which 
is significant compared to the error bars, in particular for the 
217 GHz map. To determine whether this dip is in fact responsi- 
ble for large deviation in the reconstruction, we take the A - 10^ 
best fit reconstruction and set f{k) = everywhere except for 
0.1 Mpc ' < k < 0.15 Mpc ' — the region where the large devi- 
ation is located — and calculated the corresponding C( spectrum. 
Fig. 17 plots the Cf residuals of this test feature, which show a 
large dip at around £ 1800, thus demonstrating that the dip 
in the Q residual of the data centred at£x 1800 is responsible 
for the large excursions in the primordial power spectrum recon- 
structions. The features in the temperature power spectrum, par- 
ticularly the broad dip at ^ a; 1800, cannot be explained at this 
level by any of the known systematics that have been propagated 
through the the full data analysis pipeline (Planck Collaboration 
VI, 2013). Nevertheless, we caution the reader that these (and 
other small) departures from the best fit ACDM spectrum may 
possibly be due to unknown systematic effects, or inaccurate 
propagation of known systematics, into the final power spectra. 
This possibility will be investigated intensively over the coming 
year using the data from the full mission. 



8. Parametric searches for primordial power 
spectrum features 

In this section we continue to investigate deviations of the pri- 
mordial power spectrum from a smooth, featureless function, in 
this case by testing a set of theoretically motivated models. 



8. 1. Models and priors 

We consider three models describing features in the primordial 
power spectrum, adding a global oscillation, a localized oscilla- 
tion, or a cutoff to the large scale power spectrum. 



Here we consider the following parameterization of the pri- 
mordial spectrum (referred to as the wiggles model): 



W) = mic) n + a„ sin 



win I — I + if 



(50) 



with amplitude frequency w, and phase ip to quantify the 
superimposed oscillations. The underlying smooth spectrum has 
the standard power law form 



n(fc) = A3 1- 



(51) 



The prior ranges for the wiggles model parameters are given 
in Table 10. The obvious prior for the phase (p is uniform over 
the interval (0, 2n). We choose a uniform prior on (a loga- 
rithmic prior on introduces considerable dependence of the 
resulting marginalized posteriors on the lower limit and does not 
contain the smooth spectrum as a limiting case). The sensitiv- 
ity to primordial wiggles is limited at high frequencies by the 
width of the transfer function (Hamann et al., 2008a) and at low 
frequencies by the requirement of at least one full oscillation in 
the observable part of the power spectrum. Since Planck data 
are sensitive to wavenumbers over a range of roughly six orders 
of magnitude, this condition implies u > 0.5. Here we restrict 
the analysis to co < 100 and assume a uniform prior Larger 
values of the frequency are theoretically possible, e.g., in axion 
monodromy models (Flauger et al., 2010), but the amplitude of 
the oscillations in the Cf s will be suppressed with respect to the 
primordial one. A comprehensive search for higher frequency 
oscillations is currently underway. 

8.1.2. Step-inflation model 

If the slow roll of inflation is briefly interrupted, for instance by 
a phase transition (Starobinsky, 1998; Hunt & Sarkar, 2004), a 
burst of resonant particle production (Chung et al., 2000), a sud- 
den turn in field space (Achucarro et al., 201 1) or a step in the in- 
flaton potential (Adams et al., 2001), a localized oscillatory fea- 
ture is superimposed on the scalar primordial power spectrum. 
We adopt the approximate parameterization for such a feature 
from a step in the potential, introduced by Adshead et al. (2012), 
with 



r^ik) = exp 



3 sinh(Aj)7f/.)Cd) 



where 



W'(x) 



-3 -H ^|cos2x + 1 15 



9 \ sin2x 
I 2x 



(52) 



(53) 



8.1 .1 . Wiggles model 

Due to the exponential growth of the scale factor during in- 
flation, a periodically recurring event in proper time which af- 
fects the amplitude of curvature perturbations would produce 
features that are periodic in In A:. This occurs, for instance, 
for non-Bunch-Davies initial conditions (Easther et al., 2001; 
Danielsson, 2002; Martin & Brandenberger, 2003; Bozza et al., 
2003), or, e.g., in the axion monodromy model (Silverstein & 
Westphal, 2008), as a consequence of instanton-induced correc- 
tions to the potential (Flauger et al., 2010). In these scenarios 
the primordial spectrum has an oscillation superimposed on an 
underlying smooth spectrum. 



As in the wiggles model, we choose a uniform prior on the am- 
plitude parameter yif (see Table 10). The parameter rjf deter- 
mines both the frequency of the feature and its location, which 
is required to lie in the observable range. The damping envelope 
of the feature is set by the ratio T/f/xj. We impose uniform priors 
on the logarithms of rjf and x^. 

8.1.3. Cutoff model 

A number of models were suggested to explain the apparent lack 
of power in the quadrupole and octupole of the WMAP tem- 
perature power spectrum. Typically in these models, the onset 
of a slow-roll phase coincides with the time when the largest 
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and cutoff model parameters. 



observable scales exited the horizon during inflation. This nat- 
urally suppresses the primordial power spectrum at large scales 
(see, e.g., Sinha & Souradeep, 2006). We consider a phenomeno- 
logical parameterization of a cutoff proposed in Contaldi et al. 
(2003), given by 



(54) 



We apply uniform priors on A^, which determines the steepness 
of the cutoff, and on the logarithm of the cutoff scale k^. 

8.2. Method 

To achieve the necessary numerical precision for models with 
features in the primordial spectra, we modify the standard set- 
tings of the CAMB numerical code in order to calculate C( at 
each { rather than interpolating and refine the grid in wavenum- 
ber for the numerical integration. These changes significantly 
slow down the computation. In the models considered here, 
the likelihood function has characteristics that make sampling 
difficult, such as extended plateaus and multiple isolated max- 
ima, rendering the Metropolis-Hastings algorithm inefficient. 
We therefore use the nested sampling algorithm implemented in 
the MultiNest add-on (Feroz & Hobson, 2008; Feroz et al., 
2009) to CosmoMC (Lewis & Bridle, 2002), which is also able 
to calculate the Bayesian evidence and likelihood profiles. 

The signatures of the feature models under investigation are 
rather unique and generally cannot be mimicked by other param- 
eters, which only lead to smooth variations of the power spec- 
trum (with an exception of highly tuned very low-frequency os- 
cillations that can change the acoustic peak structure). We thus 
restrict ourselves to varying only the parameters describing the 
features and keep all remaining cosmological and nuisance pa- 
rameters fixed to their ACDIVI best-fit values.''* 

8.3. Results 

For all three models we find that including these additional fea- 
tures improves the quality of the fit with respect to a pure power 
law spectrum. For the Planck+WP data, we show the best-fit pri- 
mordial curvature power spectra and temperature angular power 



An a-posteriori maximization of the likelihood in a narrow parame- 
ter range around the best-fit feature model parameters, including a varia- 
tion of all remaining cosmological and nuisance parameters, shows that 
the change in the best-fit is merely 0(1) and hence does not affect 
our conclusions. 
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Fig. 18. Top: Best fit primordial spectrum of curvature pertur- 
bations for the power law (black), wiggles (red), step-inflation 
(green), and cutoff (blue) models. Bottom: Residuals of the tem- 
perature angular power spectrum. Note that the scale of the ver- 
tical axis changes at £ = 50. Inset: Zoom on the region of the 
first acoustic peak. 



spectrum residuals in Fig. 18, and report the best fit parameter 
values in Table 1 1 . Since in all three cases the likelihood func- 
tions do not tend to zero in all directions of the respective param- 
eter spaces, the Bayesian quantities (i.e., posterior distributions 
and Bayes factors) depend considerably on the choice of prior. 
For this reason, we also quote two prior-independent quantities, 
the effective x' (i.e., -2A In = 2 In - 2 In X^;™^) and 
the profile -2A In X-max as a function of selected model param- 
eters plotted alongside the marginalized posteriors in Fig. 19, 
which illustrates the unconventional shape of the likelihood 
functions. 

For the wiggles model, oscillations around the first acous- 
tic peak and in the 700 < { < 900 range improve the fit to the 
data, whereas for the best fit step-inflation model the spectrum 
between the Sachs-Wolfe plateau and the first acoustic peak is 
fit better Quantitatively, the cutoff model improves the fit only 
modestly, with Axl^ ~ 3, but both the wiggles and step-inflation 
models lead to a larger improvement, with A^lg ~ 10, at the cost 
of three new parameters. Already for pre-Planck data, improve- 
ments of ^xlff ~ 10 have been reported in related analyses (e.g., 
Peiris et al. 2003; Martin & Ringeval 2004; Elgar0y et al. 2003; 
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Covi et al. 2006; Meerburg et al. 2012; Benetti et al. 2013; Peiris 
et al. 2013). Note that in the step-inflation model, the best-fit 
does not coincide with the maximum of the marginalized pos- 
terior probability, indicating that some degree of fine-tuning is 
necessary to reach the maximum of the likelihood. The maxi- 
mum of the marginalized posterior at In (//f /Mpc) » 7.2 actually 
reproduces the feature at /' « 20 -40 found previously in WMAP 
data (Peiris et al., 2003). It is also worth pointing out that the 
secondary peak at In (?7f/Mpc) » 4 corresponds to a feature at 
multipoles i ^ 1800, i.e., where the analysis of Sect. 7 found a 
feature. However, the model is not able to describe this feature 
very well, yielding only an improvement of ^xlff ~ 3. 

Whether or not these findings can be considered statistically 
significant or arise simply from over-fitting noisy data is not a 
trivial question (see, for instance, the discussion in Bennett et al. 
(2011)). From a frequentist statistics point of view, an answer 
would require the rather involved procedure of repeating the 
analysis on a large set of simulations. In designing the test statis- 
tic, special care would need to be taken in making sure to take 
into account the look elsewhere effect (i.e., the fact that a partic- 
ular observed anomaly may be very unlikely, whereas the prob- 
ability of observing some anomaly may be much larger). From a 
Bayesian statistics point of view, it is the Bayesian evidence that 
can tell us how probable the extended models are, compared to 
the baseline power-law primordial power spectrum. 

For the models considered here, and for our choice of 
prior probabilities, the Bayesian evidence in fact favours, albeit 
weakly, the simple power-law spectrum over the more complex 
models. The reason is that the Bayesian evidence punishes a lack 
of predictivity in these models: most of the volume of their pa- 
rameter space is not compatible with the data; a good match to 
observations is obtained only in a small region. Nonetheless, the 
observed features remain interesting, since, if they are real, they 
will also leave traces in other observables. Most notably, in the 
ii-mode polarization spectrum, where the signatures of features 
in the primordial spectrum are actually less washed out than in 
the temperature spectrum (Mortonson et al., 2009) — the forth- 
coming Planck polarization data will prove very useful in this 
regard. Additionally, since strong deviations from power law be- 
havior typically indicate non-linear physics, these models gener- 
ically also predict a non-Gaussian signal which is potentially ob- 
servable in the bispectrum (Planck Collaboration XXIV, 2013). 
However, the best-fit wiggles and step-inflation models have os- 
cillations with too high a frequency to be accessible to bispec- 
trum analysis at present. 



9. Combined analysis with P/anc/f /nl constraints 
for single field inflation 

In the previous sections we have analysed inflationary models 
with a canonical kinetic term. This led to the tensor-to-scalar 
consistency condition, requiring - -r/8. It is interesting to 
consider more general classes of inflationary models charac- 
terized by a non-standard kinetic term or more general higher 
derivative operators. An interesting subclass of these models 
are those in which the Lagrangian is a general function of the 
scalar inflaton field and its first derivative: £ - P{(p, X), where 



A more general extension is provided by the 



so-called effective field theory of inflation (Cheung et al., 2008), 
which has a richer phenomenology. 

We restrict our analysis to the first class of models (Garriga 
& Mukhanov, 1999; Chen et al., 2007), which includes k- 
inflation models (Armendariz-Picon et al. 1999; Garriga & 
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Table 11. Improvement in fit and logarithm of the Bayes factor 
with respect to power law ACDM and best fit parameter val- 
ues for the wiggles, step-inflation, and cutoff models. The larger 
In Box, the greater the preference for a featureless power law 
spectrum. 



Mukhanov 1999), and Dirac-Born-Infield (DBI) models intro- 
duced in the context of brane inflation (Silverstein & Tong 2004; 
Alishahiha et al. 2004). In this class of models inflation can take 
place with a steep potential, or it can be driven by the kinetic 
term. One of the main features of inflationary models with a 
non-standard kinetic term is that the inflaton fluctuations can 
propagate at a sound speed Cj < 1 . As shown in previous anal- 
yses (e.g., Peiris et al., 2007; Powell et al., 2009; Lorenz et al., 
2008; Agarwal & Bean, 2009) there are strong degeneracies be- 
tween the parameters determining the observable power spec- 
tra. Constraints on primordial non-Gaussianity can help break 
this degeneracy, and we show how Planck's combined measure- 
ment of the power spectrum and the non-linearity parameter /nl 
(Planck Collaboration XXIV, 2013) improves constraints on this 
class of models. 

In models with a non-standard kinetic term the sound speed 
of the inflaton is given by cl - P^x/{P,x + 2Xf xx) (Garriga 
& Mukhanov, 1999), so that in the canonical models, where 
P((f>, X) = V(<p) - X, one finds Cs = 1 , while in general a non- 
trivial Cs < 1 corresponds to deviations from this standard case. 
Therefore, in these models, new parameters — like the sound 
speed and its running — enter the expressions of the inflationary 
observables. For the running of the sound speed it is useful to 
define an additional slow-roll parameter 



s - 



c,H 



(55) 



For values of the slow-roll parameters much less than unity, 
the leading order scalar power spectrum is modified (Garriga & 
Mukhanov, 1999) to 



1 //^ 
pi 



(56) 



which is evaluated at kc^ - aH. The scalar spectral index gets an 
additional contribution from the running of the sound speed. 



- 1 = -2ei - €2- s . 



(57) 



The gravitational sector remains unaltered by the non-trivial 
inflaton sound speed, retaining the same form as for the stan- 
dard slow-roll models. Therefore the usual consistency relation 
is modified to r ^ -8ntCs with «t = -2e\ as usual (Garriga & 
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Fig. 19. Marginalized posterior probability {red) and profile 
-2A In -Cmax (black) for selected parameters of the wiggles {top), 
step-inflation {middle), and cutoff model {bottom). 

Mukhanov, 1999). A more accurate relation that we employ in 
this analysis is given by 

r = 16eic,'"'' . (58) 

This accounts for the difference in freeze-out between the scalar 
and tensor perturbations (Peiris et al., 2007; Powell et al., 2009; 
Lorenz et al., 2008; Agarwal & Bean, 2009), taking place at 



sound-horizon crossing (at kc^ — aH) for the scalar fluctuations 
and at horizon-crossing {k = aH) for the tensor modes. 

Limiting ourselves to the predictions at lowest-order in the 
slow-roll parameters, there are clearly degeneracies between the 
parameters {As,c^,ei,e2, s), which make the constraints on the 
inflationary power spectra observables less stringent in terms 
of these microscopic parameters. However, for models where 
the inflaton field has a non-standard kinetic term with Cj <s; 1, 
a high level of primordial non-Gaussianity of the scalar per- 
turbations is generated (see, e.g., Chen et al., 2007). In these 
models primordial non-Gaussianity is produced by the higher 
derivative interaction terms that arise when expanding the ki- 
netic part of the Lagrangian, P{(p, X). The amplitude of the non- 
Gaussianity, defined by the nonlinearity parameter /nl, receives 
two dominant contributions, arising from the inflaton interac- 
tion terms {S(p)^ and 6(f){'V64>)^. Each of them produces non- 
Gaussianity shapes similar to the the so-called "equilateral" 
type (Babich et al. 2004), i.e., a signal that peaks for equilateral 
triangles k[ - k2 - k^,. However, they are sufficiently different 
that the total signal can be very different from the equilateral 
one (Senatore et al., 2010). The non-linearity parameter of the 
second interaction term is /nl = (85/324)(l - c^"^), while the 
other is determined by a second independent amplitude (Chen 
et al., 2007; Senatore et al., 2010). Constraints on the primor- 
dial non-Gaussianity, presented in the companion paper Planck 
Collaboration XXIV (2013), thus allow us to construct a lower 
limit of the sound speed c^. This helps reduce degeneracies in 
the parameter space of inflationary models with non-standard ki- 
netic terms. In particular, without the limits on the sound speed 
coming from the constraints on primordial non-Gaussianity, it 
is not possible to derive an upper limit on the parameter ei , be- 
cause the relation between the tensor-to-scalar ratio and ej also 
involves the sound speed (see, e.g., Eq. 58). 

In this paper, we consider three cases. One is a general anal- 
ysis as described above, where we focus on the simplest case of 
a constant speed of sound with s - Q. From the Planck limits 
on primordial non-Gaussianity in general single field models of 
inflation (Planck Collaboration XXIV, 2013), the most conser- 
vative constraint on the sound speed is 

Cs > 0.02 (95% CL) . (59) 

In this large parameter space, we assume a uniform prior 
0.02 < Cs < 1 in Eq. 58 within the HFF formalism described 
in the Appendix. We show the joint constraints on ei and €2 in 
Fig. 20. By including the 95% CL constraint on from Eq. 59, 
Planck+WP constrain ei < 0.053. Such constraints can be com- 
pared with the restricted case of c^ - 1, also shown in Fig. 20, 
with ei < 0.008 at 95% CL. 

The other two cases we analyse correspond to two specific 
models where the inflaton has a non-standard kinetic term. The 
degeneracy between the different slow-roll parameters is broken 
because these models specifically predict that s = 0, or ^ oc ^2. 
As an example, we first consider the case where the action takes 
the Dirac-Born-Infield (DBI) form 

p{<p, X) = -/(0)-i Vi - VWx + m-' - y(^) . (60) 

Here y(0) is the potential and f{(p) is the warp factor deter- 
mined by the geometry of the extra dimensions. For DBI mod- 
els a stronger bound on Cs is derived (Planck Collaboration 
XXIV, 2013): c, > 0.07 at 95% CL. With the uniform prior 
0.07 < Cs < 1 and 5 = 0, Planck + WP constrain ei < 0.042 at 
95% CL. 
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Fig. 20. Marginalized joint 68% and 95% CL for (ei ,62) for 
Planck + WP data, comparing the canonical Lagrangian case 
with Cs = 1, to the case with varying with a uniform prior 
0.02 < Cs < 1, derived from the Planck non-Gaussianity mea- 
surements. 

An important case is f(^) ^ A/^'^ (for details, see Silverstein 
& Tong 2004; Alishahiha et al. 2004; Chen et al. 2007, and ref- 
erences therein). There are two possibilities. First, in ultraviolet 
(UV) DBI models, the inflaton field moves under a quadratic 
potential V{^) ^ ifp-(p^/2 from the UV side of the warped space 
to the infrared side, with m » M^i/ y/A. It is known that this 
case is already at odds with observations, if theoretical internal 
consistency of the model and constraints on power spectra and 
primordial non-Gaussianity are taken into account (Baumann & 
McAllister, 2007; Lidsey & Huston, 2007; Bean et al., 2007; 
Peiris et al., 2007; Bean et al., 2008). It is therefore interesting 
to look at the other case, namely infrared (IR) DBI models (Chen 
2005b, a) where the inflaton field moves from the IR to the UV 
side, and the inflaton potential is 

V(0) = Vo - , (61) 

with a wide range of values allowed for p in principle, 
0.1 < /3 < 10^ (Bean et al., 2008). Here, we focus on a mini- 
mal version of the IR DBI models, where string effects are ne- 
glected, so that the usual field theory computation of the pri- 
mordial curvature perturbation holds. For IR DBI models ac- 
counting for such effects and a more involved treatment of the 
dynamics, see Chen (2005a,c); Bean et al. (2008). In this mini- 
mal IR DBI model, one finds (Chen, 2005c; Chen et al., 2007) 
c, =^ (J3NJ3)-\ ns - 1 = -4/A^„ dnjdlnk = (in this 

model one can verify that s ^ ^ ^2/3). Here, primordial 

non-Gaussianity of the equilateral type is generated with an am- 
phmde /D^^^ = -(35/108) [(/3^N^/9) - 1]. 

If we consider 60 < Nt < 90, then the predicted spectral 
index lies within the range 0.93 < n^ < 0.96, which is con- 
sistent with the Planck measurement of the spectral index at 
the 3cr level, for > 60. The constraints on non-Gaussianity 
give /™' = 11 ± 69 at 68% CL (Planck Collaboration XXIV, 
2013). Combining these constraints with the power spectrum 
constraints, marginalizing over 60 < Nt < 90, we obtain 

/3<0J (95% CL), (62) 

This strongly limits the allowed parameter space of these mod- 
els. 



As a final example, we consider a class of power-law k- 
inflation models characterized by the Lagrangian (Armendariz- 
Picon et al., 1999) 

P((I>,X) = Ii::^l(-X + X2). (63) 

In this case, for small values of y one finds: ^ y/8, "Pr = 
IH^/iSjcSn^M^^Xk/ko)-^^, n.s - 1 = -3r- The sound speed 
is a constant (s - 0), with constant y. The primordial non- 
Gaussianity in this model has an amplitude /^l" ' ~ -170/(81y). 
Therefore, all the inflationary observables depend essentially on 
a single parameter, y. Imposing a prior of < y < 2/3, from the 

non-Gaussianity constraint /j^l ~ - ^^"^ (Planck 
Coflaboration XXIV, 2013) we obtain y > 0.05 at 95% CL. At 
the same time, our measurement of the spectral index constrains 
0.01 < y < 0.02 at 95% CL. This class of /t-inflation models 
is therefore excluded by the combined constraints on primordial 
non-Gaussianity and the power spectrum. 



10. Isocurvature modes 

10.1. Theoretical background 

In this section we explore the constraints imposed by Planck on 
scenarios where the primordial cosmological perturbations were 
not entirely adiabatic. These scenarios also include isocurvature 
modes, possibly correlated among themselves as well as with 
the adiabatic mode. The adiabatic mode is characterized by the 
property that at very early times the Universe obeyed a com- 
mon, spatially-uniform equation of state, and all components ini- 
tially shared a common velocity field. For the adiabatic mode the 
density perturbations in the various components (e.g., baryons, 
CDM, photons, neutrinos) are locked together. Here baryons in- 
clude their accompanying leptons, assumed tightly coupled to 
maintain charge neutrality. 

Isocurvature modes arise from spatial variations in the equa- 
tion of state or from relative velocities between the compo- 
nents. To analyse how the CMB perturbations were imprinted, 
it is most convenient to define isocurvature modes at a suffi- 
ciently late time such that the relevant components, according to 
our present best understanding, consisted of baryons, photons, 
CDM, and neutrinos. Under this hypothesis, in addition to the 
adiabatic mode there are four possible non-decaying isocurva- 
ture modes: the baryon, CDM, and neutrino density isocurvature 
modes, and the neutrino velocity isocurvature mode (see, e.g., 
Bucher et al. (2000) for a discussion and further references). 

The impact of isocurvature modes on the CMB was first 
studied in detail by Peebles & Yu (1970) and Efstathiou & Bond 
(1986, 1987), who contemplated the possibility that isocurvature 
perturbations rather than adiabatic perturbations were the sole 
source of cosmological perturbations. Linde (1985), Polarski & 
Starobinsky (1994), Linde & Mukhanov (1997); and Garcfa- 
Bellido (& Wands (1996) pointed out various scenarios in which 
isocurvature perturbations could be generated within the con- 
text of inflation. Bucher et al. (2000) carried out a systematic 
study of isocurvature modes from a phenomenological perspec- 
tive, pointing out the relevance of two additional modes (the neu- 
trino density and velocity modes). Lyth & Wands (2002), Moroi 
& Takahashi (2001) and Bartolo & Liddle (2002) studied an in- 
teresting so-called curvaton scenario, in which adiabatic fluctu- 
ations from inflation contribute negligibly but quantum fluctua- 
tions in a transverse direction modulate the density of decaying 
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particles, leading to isocurvature perturbations correlated with 
the adiabatic mode. 

Several authors have studied the constraints on isocurva- 
ture modes imposed by previous microwave background exper- 
iments, including Stompor et al. (1996), Langlois & Riazuelo 
(2000), Amendola et al. (2002), Peiris et al. (2003), VaUviita & 
Muhonen (2003), Bucher et al. (2004), Moodley et al. (2004), 
Beltran et al. (2004), Kurki-Suonio et al. (2005), Dunkley et al. 
(2005), Bean et al. (2006), Trotta (2007), Keskitalo et al. (2007), 
and Komatsu et al. (2009). A more complete set of references 
may be found in Valiviita et al. (2012). 

Before proceeding we must define precisely how to char- 
acterize these isocurvature modes on superhorizon scales dur- 
ing the post-entropy-generation epoch, during which we assume 
that the stress-energy content of the Universe can be modelled as 
a multi-component fluid composed of baryons, CDM particles, 
photons, and neutrinos. If we assume that the evolution of the 
Universe during this epoch was adiabatic (used here in the sense 
of thermodynamically reversible), then the entropy per unit co- 
moving volume is conserved and serves as a useful reference 
with respect to which the abundances of the other components 
can be expressed. 

The baryon isocurvature mode may be expressed in terms of 
fractional fluctuations in the baryon-to-entropy ratio, which is 
conserved on superhorizon scales during this epoch. The CDM 
and neutrino density isocurvature (NDI) modes may be defined 
analogously. The neutrino velocity isocurvature (NVI) modes 
refer to fluctuations in the neutrino velocity relatively to the 
average bulk velocity of the cosmic fluid. For the CMB, the 
baryon and CDM isocurvature modes yield almost identical an- 
gular spectra, because the deficit of one is balanced by an ex- 
cess of the other, so we do not consider them separately here. 
In this way the primordial isocurvature modes may be defined 
as dimensionless stochastic variables Jcdi, -^ndi, -^nvi, like the 
variable K describing the adiabatic mode.'"" In this basis, the 
CDI mode can be seen as an effective isocurvature mode, en- 
coding both CDM and baryon isocurvature fluctuations through 
JcDi + (i^b/i^c)i'Bi (Gordon & Lewis, 2003). 
Within this framework, Gaussian fluctuations for the most 
general cosmological perturbation are described by a 4 x 4 pos- 
itive definite matrix-valued power spectrum of the form 
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(64) 



Following the conventions used in CAMB (Lewis & Bridle, 2002; 
Lewis, 201 1) and CLASS (Lesgourgues, 201 1; Bias et al., 201 1), 
the primordial isocurvature modes are normalized as follows 
in the synchronous gauge: for the CDI mode, Vnik) is the 
primordial power spectrum of the density contrast difference 
^Pcdm/pcdm; for the NDI mode it is that of Spvlpv', and for the 
NVI mode, that of the neutrino velocity Vy times 4/3"^. 

If isocurvature modes are present, the most plausible 
mechanism for exciting them involves inflation with a multi- 
component inflaton field. To have an interesting spectrum on 
the large scales probed by the CMB, isocurvature modes require 

The symbol S is sometimes used in the literature to denote the 
isocurvature modes, also known as entropy perturbations. To prevent 
confusion we avoid this terminology because isocurvature modes are 
unrelated to any notion of thermodynamic entropy. 

or in other words, of the neutrino perturbation dipole, 
F^i = 4e„/(3yt) in the notation of Ma & Bertschinger (1995). 





Fig. 21. Ctt anisotropy shape for the three isocurvature modes. 
Top: The shapes of the CDM isocurvature mode, neutrino den- 
sity isocurvature mode and neutrino velocity isocurvature mode 
are shown together with the adiabatic mode. The modes have 
the same amplitude parameters {Vkr for the adiabatic mode and 
Vii for each isocurvature mode). Bottom: the narrower multi- 
pole range illustrates the relative phases of the acoustic oscilla- 
tions for these modes. 



long-range correlations, and at present inflation with a multi- 
component inflaton provides a well-motivated scenario for es- 
tablishing such correlations. Inflation with a single component 
scalar field can excite only the adiabatic mode. In models of 
inflation with light (compared to the Hubble expansion rate) 
transverse directions, the scalar field along these transverse di- 
rections becomes disordered in a way described by an approxi- 
mately scale-invariant spectrum. If the inflaton has light com- 
ponents, there are (A^- 1) potential isocurvature modes during in- 
flation. Whether or not the fluctuations along these transverse di- 
rections are subsequently transformed into the late-time isocur- 
vature modes described above depends on the details of what 
happens after inflation, as described more formally below. 

As explained for example in Gordon et al. (2001), 
Groot Nibbelink & van Tent (2000, 2002) and Byrnes & Wands 
(2006), for inflationary models where the inflaton follows a 
curved trajectory, correlations are generically established be- 
tween the isocurvature and curvature degrees of freedom. To 
lowest order in the slow-roll approximation, this leads to a situa- 
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tion where the adiabatic perturbation is the sum of several com- 
ponents each of differing spectral index. 

The post-inflationary evolution determines how the isocur- 
vature fluctuations generated during inflation transmute into the 
three specific isocurvature modes studied here. Little is known 
about the details of what happens during the epoch of entropy 
generation, but to linear order we may express how the fields 
■Rinf (i.e., the curvature perturbation at the end of inflation) and 
the transverse components of the inflaton field cri, . . ., cr^-i (i.e., 
the components orthogonal to the slow-roll direction) transform 
into curvature perturbations and late-time isocurvature modes at 
the end of the epoch of entropy generation as the following linear 
transformation: 



^out 



1 

Mm 



O-A 



(65) 



where fl = 1,...,4 = BI, CDI, NDI, NVIandA = 1, . . . , (A? - 1) 
labels the transverse components of the component inflaton 
field. Physically, the fluctuations along the transverse directions 
modulate particle production during the epoch of entropy gener- 
ation. 

The neutrino density isocurvature can be excited in much the 
same way as the CDM and baryon isocurvature mode because 
at least within the standard electroweak model, in which there 
are no leptonic flavour changing processes, L^, L^, and Lj are 
separately conserved. Known non-perturbative processes such 
as the sphaleron can turn trade lepton and baryon asymmetries 
with each other and alter flavour asymmetries, but they cannot 
erase such asymmetries altogether Plausible models generating 
the neutrino density mode are therefore possible (Bucher et al., 
2000; Gordon & Malik, 2004), but as for the neutrino velocity 
isocurvature mode, to date no plausible generation mechanism 
has been put forth. 

This extension of the adiabatic ACDM model to non- 
adiabatic initial conditions represents an important test of in- 
flation. Single field inflation can produce only adiabatic pertur- 
bations, since exciting isocurvature perturbations requires addi- 
tional degrees of freedom during inflation. Therefore a detection 
of primordial isocurvature perturbations would point to more 
complicated models of inflation. 

10.2. Adiabatic witli one isocurvature mode and free 
spectral indices 

In this paper we investigate three of the four possible isocurva- 
ture modes of the ACDM scenario, since the baryon and CDM 
isocurvature perturbations are indistinguishable in the CMB an- 
gular power spectra. The CDM, neutrino density, and neutrino 
velocity isocurvature perturbations lead to different power spec- 
tra for CMB anisotropics, as shown in Fig. 21. We limit our- 
selves to studying one isocurvature mode at a time, in the pres- 
ence of a curvature perturbation. More general combinations 
with two or three isocurvature modes may be contemplated, but 
without the Planck high frequency polarization likelihood, it is 
difficult to constrain this scenario, so we postpone a discussion 
of this case to the next release. 

Theoretically, one expects the power spectra of the isocurva- 
ture modes and their correlations to exhibit near but not neces- 
sarily exact scale invariance. As a general test of adiabaticity, it 
is nevertheless interesting to compare a more general model to 
the Planck data, assuming that the adiabatic, isocurvature, and 
cross-correlation spectra obey a power law with free spectral 
indices. Blue values of the spectral indices are particularly in- 
teresting from the point of view of testing adiabaticity because 



the acoustic peaks arising from two of the isocurvature modes 
are out of phase with the adiabatic peaks by roughly (n/T) near 
the first acoustic peak. This is not true for the neutrino velocity 
isocurvature mode however. 

In the literature, models with one isocurvature and the adi- 
abatic mode (possibly correlated) are often parameterized by 
specifying the 2x2 correlation at a certain pivot scale ko through 
the components ^Prj, 'Pii along with their respective spec- 
tral indices n^n, nm, nil (e.g., Amendola et al., 2002; Beltran 
et al., 2004). We do not follow this approach here because in 
the absence of a statistically significant detection, the posterior 
distributions for the spectral indices are difficult to interpret and 
sensitive to how the prior is chosen. We instead adopt a parame- 
terization where !Pab is specified at two scales k -k\ and k -ki 
and interpolated geometrically according to the relation'^ 



!Pab(^)=exp 



In(fc) - ln(fc2) 
In(^i) - ln(^2) 



^ / In(fc) 



In(fci) 



\Mki)-Mkx) 



In(^^fb) 



(66) 



where a,b - I and I -Ica\, -^ndi> or Jnvi- We set 
fci = 2 X 10"^^ Mpc"' and A:2 = 0.1 Mpc"', so that \kx,kT\ spans 
most of the range in k constrained by Planck data. A uniform 
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prior for the components V, 

assumed, where auto-correlation amplitudes !P^^, ^jj^ 
!Pjj are positive, although the cross-correlation amplitudes V. 
'P'^i, may take both signs subject to the constraints 
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«j' 



m' II' 



KR' II' 



(67) 



to ensure positive definiteness. For the logarithm of the off- 
diagonal elements in Eq. 66 to be real, we must have 
'^Kr'Pm > 0. This Ansatz can be trivially modified to admit 
the case < by inserting appropriate minus signs, 

but this parameterization does not allow the case where the sign 
of the correlation changes. In practice we deal with this by as- 
suming a uniform prior not on !P^], but on its absolute value. 



3(2) 



sign(!P^']) X The constraints 



and then we impose V^^j 

in Eq. 67 ensure that det (!Pab(^)) is positive definite within the 
interval {ky,k2\, but generically positive definiteness is violated 
sufficiently far outside this interval, either for very small or very 
large k. Where this happens we reduce the magnitude of Vkj so 
that there is either total correlation or anti-correlation. The kinks 
thus introduced lie outside the range [^i , ^2] ■ Within this range, 
the spectral indices n-R^, n^j, njj are scale-independent. Finally, 
our sign conventions are such that positive values for cor- 
respond to a positive contribution of the cross-correlation term to 
the Sachs-Wolfe component of the total temperature spectrum. 

When the constraining power of the data is weak, a crucial 
question is to what extent the posterior distribution results from 
the data rather than from the prior distribution. The parameter- 
ization above is not the only one that could have been adopted, 
and other possible priors are typically related by a non-constant 
Jacobian. For each model, we indicate the log-likelihood for the 
best fit model, in order to allow model comparison. 

The Planck+WP results for the three isocurvature modes us- 
ing this two-scale parameterization are shown in Fig. 22 and 



" Although the models spanned by the one-scale and two-scale pa- 
rameterizations are the same, the priors for these parameterizations are 
related by a non-constant Jacobian and therefore do not coincide. 
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included in the summary Table 12. The power spectra fmik), 
V-Riik), and Vnik) are normalized according to the primordial 
values of the fields 'R(x) and J(x) defined above. It is interest- 
ing to consider how much isocurvature is allowed expressed as a 
fraction of the power in three bands spanning the CMB temper- 
ature spectrum observed by Planck. To this end, we define the 
following derived quantities 



£^KK(^min, ^max) 



ffKj(^min, ^max) - 



(^^)^^(^min? ^max) 
(^-^)tot(^min' ^max) 
(^^)jj(^min? ^max) 
(^^)tot(^min' ^max) 
(^^)^j(^min5 ^max) 
(^^)tot(^iiiin5 ^max) 



(68) 
(69) 
(70) 



where 



{^T)li(^^, {^) = J] (2^ + \)Ciy (71) 

The 95% confidence limits from the one-dimensional posterior 
distributions for these fractional contributions in the full range 
(^min,^max) = (2,2500) are shown in Table 12. The range of al- 
lowed values for okkC^-, 2500) is a measure of the adiabaticity of 
fluctuations in the CMB. The posterior distributions of the frac- 
tions ajj, OKI in three multipole ranges are shown in Fig. 23. 
We also report the primordial isocurvature fraction, defined as 



PUk) 



riiik) 



'Pm(k)+'Pii(k) 



(72) 



at three values of k. Table 12 also shows the effective - 
-21nXmax for all models, compared to the minimal six- 
parameter ACDM model. In Fig. 24 we show the ratio of temper- 
ature spectra for the best-fit mixed model to the adiabatic model. 
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Fig. 22. Two dimensional distributions for power in isocurvature 
modes, using Planck+WP data. 
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Fig. 23. Fractional contribution of isocurvature modes to the 
spectrum. We show the distributions ci'jj(2, 20), a'Kj(2,20), 
Of jj(2 1,200), a?jj(2 1,200), a7j(201,2500), aKj(201,2500) de- 
fined in Eq. 70 for the CDI, NDI, NVI modes, constrained with 
Planck+WP data. 




Fig. 24. Temperature spectrum of best-fit models with a mix- 
ture of adiabatic and isocurvature modes. Top: spectrum of the 
best-fit mixed models relative to that of the pure adiabatic model. 
Bottom: zoom on the Sachs-Wolfe plateau of the best-fit temper- 
ature spectrum D( - {[{{ + 1)/27t]CJ^, for each of the three cases 
plus the pure adiabatic model, shown together with Planck low-{ 
data points. 



The results for a'????(2, 2500) show that the nonadiabatic con- 
tribution to the temperature variance can be as large as 7% (9%, 
5%) in the CDI (NDI, NVI) model (95% CL). These results are 
driven by the fact that on large scales, for I < 40, the Planck 
data points on average have a slightly smaller amplitude than 
the best-fitting ACDM model. Hence the data prefer a signifi- 
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Uncorrelated, njj = 1, ("axion") 
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Table 12. Isocurvature mode constraints. For each model, we report the 95% CL upper bound on the fractional primordial contribu- 
tion of isocurvature modes at three co-moving wavenumbers (A:iow = 0.002 Mpc"', ^mid = 0.05 Mpc"', and A:high - 0.10 Mpc"'), and 
the 95% CL bounds on the fractional contribution ajj and a^j^j to the total CMB temperature anisotropy in the range 2 < £ < 2500. 
We also report -2AlnXmax for the best-fitting model in each case, relative to the best-fit 6-parameter ACDM model, with the 
number of additional parameters An. In the Gaussian approximation, -2AlnXmax corresponds to h.x^- The general models have six 
parameters that specify the primordial correlation matrix at two scales k\ and ki, thus allowing all spectral indices to vary (so, four 
parameters more than the pure adiabatic model). 



cant amount of anticorrelated isocurvature modes, leading to a 
reduction of amplitude of the Sachs-Wolfe plateau and to a de- 
crease of the effective by up to 4.6'**. This situation explains 
the rather loose bounds on the derived parameter aii{2, 20), as 
shown in Fig. 23. 

A comparison of and P^jj shows that best-fitting models 
have an isocurvature spectral index n jj close to 1 .7 for CDI, 1 . 1 
for NDI and 1 .0 for NVI modes. 

For CDI and NDI, the amplitude of acoustic peaks quickly 
decreases with increasing {, so the constraints are entirely driven 
by small ^'s. Since the same value of the primordial amplitude 
!P*?j leads to different plateau amplitudes for the two isocurva- 
ture models (see Fig. 21), the bounds on !Pjj and are con- 
sistently stronger for CDI than for NDI. For NVI, the acous- 
tic peak amplitude is larger than the plateau amplitude. In NVI 
models, the data cannot allow for a too large amplitude of corre- 
lated isocurvature modes at small C, because the total spectrum 
would be distorted at larger t. This possibility is strongly dis- 
favoured by the data, which is consistent with the peak location 
predicted by a pure adiabatic model. Hence in the NVI case we 
obtain slightly stronger bounds and a smaller reduction of the 
effective 

The fact that the data prefer models with a significant con- 
tribution from CDI or NDI modes should be interpreted with 
care. The detection of a shift in the phase of acoustic oscilla- 
tions would bring unambigous evidence in favour of isocurva- 
ture modes. With Planck data, we are not in this situation. The 
evidence is driven by a small deficit of amplitude in the Sachs- 
Wolfe plateau, that could have several different possible explana- 
tions (such as a deficit in the large-scale primordial power spec- 
trum, as already seen in the previous sections). However, multi- 
field inflationary scenarios can produce the mixture of curvature 
and isocurvature fluctuations which we have found to provide a 
good fit to the Planck data. 



For the three general models, the posterior distribution is actually 
multimodal. Here we are referring to models contributing to the main 
peak in the posterior, with the highest maximum likelihood. There is 
another peak with a smaller maximum likelihood, appearing in Fig. 23 
as a small bump for positive values of the cross-correlation amplitude. 
In this paper, we do not carry out a separate investigation for models 
contributing to this secondary peak. 



10.3. Special cases 

The six-parameter models of the previous subsection, includ- 
ing one isocurvature mode and the adiabatic mode, make no as- 
sumptions about the spectral indices of each mode, or the degree 
of correlation between the isocurvature mode and the adiabatic 
mode. This leads to a large number of additional degrees of free- 
dom. There are either theoretical and phenomenological motiva- 
tions for limiting the values of the parameters to certain values, 
leading to specific cases with just one more degree of freedom 
with respect to the adiabatic case. The results are reported in 
Table 12, for uncorrelated perturbations with njj - 1, and fully 
correlated or anti-correlated perturbations with njj - n^n- As 
for the general case, anti-correlated isocurvature perturbations 
slightly improve the fit to Planck data. In the following we con- 
sider the implications of our results for two important scenarios, 
the axion and curvaton scenarios. 



10.3.1. Constraints on axion isocurvature 

The axion field was proposed to solve the strong CP problem 
and constitutes a well-motivated dark matter candidate (see, 
e.g., Preskill et al. (1983), Turner (1990), Peccei (2008), Sikivie 
(2008), Raffelt (2008), and Kim & Carosi (2010)). The axion is 
the Goldstone boson of the broken Peccei-Quinn (PQ) symme- 
try. Under certain assumptions, the axion field may induce sig- 
nificant isocurvature perturbations (Turner et al., 1983; Axenides 
et al., 1983; Steinhardt & Turner, 1983; Linde, 1984, 1985; 
Seckel & Turner, 1985; Kofman, 1986; Lyth, 1990; Linde & 
Lyth, 1990; Turner & Wilczek, 1991; Linde, 1991; Lyth, 1992). 
If inflation takes place after PQ symmetry breaking, the quantum 
fluctuations of the inflaton are responsible for primordial curva- 
ture perturbations, while those of the axion field generate pri- 
mordial entropy perturbations. After the QCD transition, when 
one of the vacua becomes preferred giving the axion field a mass, 
the axions behave as cold dark matter. This way of producing ax- 
ionic dark matter is called the misalignment angle mechanism. 
In such a scenario, the CMB anisotropy may include signifi- 
cant power from CDM isocurvature fluctuations. In that case, the 
fraction ySiso = 'Pii/CPm + 'Pii) of CDM isocurvature modes is 
related to the energy scale of inflation, Hi„f, through (Lyth, 1990; 
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Beltran et al., 2007; Bae et al., 2008; Hamann et al., 2009) 



^^inf - 



0.96 X 10^ GeV / /Jj^o \'^^ / ^ 



1 0.04/ \ 0.120 



1/2 



/a 



10" GeV 



0.408 



(73) 

where Q.^ is the relic axion density, the fraction of CDM con- 
sisting of axions, and the PQ symmetry breaking scale. In this 
model, CDM isocurvature perturbations should be totally uncor- 
related with adiabatic perturbations, and have a spectral index 
njj very close to one, since in the first-order slow-roll approx- 
imation the index reads (1 - ley)- Since the sensitivity of the 
data to «jj is very limited (Beltran et al., 2007), we assume for 
simplicity that rijj = 1 . 

Within the general parametrization presented in Eq. 66, we 



can select the axion case by imposing P^'j 
condition 



<p(2) _ rpy 
^11 ~ ^11' 



3(1) 



- 0, as well as the 
(74) 



corresponding to rijj - 1. We therefore have three indepen- 



dent parameters, f^^^K "^ir sample these parame- 

ters with uniform prior distributions. The fraction j6iso(^«) with 
- 0.05 Mpc"' is then a derived parameter Since the data con- 
strain ySiso <s: 1, the relation between /3iso and fjj is nearly linear, 
so the primordial isocurvature fraction is sampled with a close- 
to-uniform prior. 

Constraints on this model are shown in Table 12. We find 



3(1) 



ySiso < 0.039 (95% CL, Planck+WP), 



(75) 



at the scale k^a = 0.05 Mpc" , with a best-fit value of zero. 
Hence there is no evidence for axion-generated isocurvature 
perturbations. This limit significantly improves the latest CMB 
bounds: at the scale k = 0.002 Mpc"', our result reads /3iso < 
0.036, to be compared with jSiso < 0.15 for WMAP 9-year alone, 
or /3iso < 0.061 for WMAP+ACT+SPT at 95% CL (Hinshaw 
et al., 2012b). This bound can be used to exclude regions in 
the parameter space composed of f^, R^, and the energy scale 
of inflation, but cannot be used to obtain a model-independent 
bound on f^. However, if we assume (i) that the PQ symmetry 
is broken during inflation, (ii) that it is not restored by the quan- 
tum fluctuations of the inflaton (which imposes //inf/(27r) < /a), 
nor by thermal fluctuations in case of a very efficient reheating 
stage, and (iii) that all of CDM consists in axions produced by 
the misalignment angle, then we can derive an upper bound on 
the energy scale of inflation as 



i/irf < 0.87 X 10^ GeV 



/a 



1011 GeV 



0.408 



(95% CL) . (76) 



10.3.2. Constraints on the curvaton scenario 

In the simplest one-field inflationary models curvature perturba- 
tions arise from quantum fluctuations in the inflaton field, but 
this is not the only way to generate them. Entropy perturbations 
may seed curvature perturbations outside the horizon (Polarski 
& Starobinsky, 1994; Gordon et al., 2001), so it is possible that 
a significant component of the observed adiabatic mode could 
be strongly correlated with an isocurvature mode. This happens 
for instance in the curvaton scenario (Mollerach, 1990; Enqvist 
& Sloth, 2002; Moroi & Takahashi, 2001; Lyth & Wands, 2002; 
Lyth et al., 2003; Gordon & Lewis, 2003). The curvaton is an 
extra light scalar field acquiring a spectrum of fluctuations on 
cosmological scales during inflation. Depending on its density 
evolution and decay history, this field could be responsible for 



part of the observed adiabatic perturbations, or all of them, or 
for a mixture of correlated adiabatic and isocurvature perturba- 
tions. 

We focus here on the simplest viable version of this sce- 
nario, in which the curvaton decays into CDM particles while 
contributing to a non-negligible fraction 



r-D 



3pci 



+ 4pn 



(77) 



of the total energy density of the Universe. If the curvaton dom- 
inates at decay time (fo - 1), its primordial fluctuations seed 
curvature perturbations equivalent to a pure adiabatic mode. If 
r/5 < 1, curvaton fluctuations are only partially converted into 
adiabatic perturbations, while CDM particles carry CDI pertur- 
bations, that are fully correlated with the adiabatic ones since 
they share a common origin. We recall that with our conventions, 
"fully correlated" means that the cross-correlation term brings a 
positive contribution to the Sachs- Wolfe component of the total 
temperature spectrum; some authors define the correlation with 
the opposite sign and and call this case "fully anti-correlated" 
(e.g., Komatsu et al., 201 1; Hinshaw et al., 2012b). In this model, 
the CDI fraction is related to ro by (Gordon & Lewis, 2003) 



IcDi _ 3(1 - rp) 

In our notation this is equivalent to 

9(l-ro)2 



^iso ~ 



r2 +9(1- rz,)2 



(78) 



(79) 



Within the general parametrization presented in Eq. 66, we can 
satisfy this case by imposing 



,(i) 

KI 
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(80) 



together with the condition 
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II 



•KR' II 



ff)(i) ' 



(81) 



corresponding to rijj - nnn. Like in the axion case, this results 



in three independent parameters "PI^^K which we sample 



3(1) 



II 



with uniform priors. The constraints for this model are shown in 
Table 12. The best-fit model is still the pure adiabatic case, and 
the upper bound 



jSiso < 0.0025 (95% CL, Planck+WP) 



(82) 



is scale independent, since the adiabatic and isocurvature tilts 
are assumed to be equal. This is a significant improvement over 
the WMAP 9-year bounds, fSi^o < 0.012 for WMAP alone, or 
Aso < 0.0076 for WMAP+ACT+SPT at 95% CL (Hinshaw et al., 
2012b). We conclude that in this scenario, the curvaton should 
decay when it dominates the energy density of the Universe, with 
ro > 0.983. 

The non-linearity parameter in the curvaton model studied 
here is (Bartolo et al., 2004c,b) 



^local 

/nl 



5 _ Std 
3 ~6~' 



(83) 



assuming a quadratic potential for the curvaton field (Sasaki 
et al., 2006). In the pure adiabatic case (ro = 1) this leads to 
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/nl ~ ~^ 1^- constraint 0.98 < < 1 then corresponds 
to -1.25 < /^l"' < -1.21. Taking into account the Planck result 
= 2.7 + 5.8 (Planck Collaboration XXIV, 2013), we con- 
clude that the Planck data are consistent with the scenario where 
the curvaton decays into CDM when it dominates the energy 
density of the Universe, and its fluctuations are almost entirely 
converted into adiabatic ones. 



11. Conclusions 

This paper establishes the status of cosmic inflation in the con- 
text of the first release of the Planck cosmological results, which 
includes the temperature data from the first 2.6 sky surveys. 
CMB polarization as measured by Planck will be the subject 
of a future release. We find that standard slow-roll single field 
inflation is compatible with the Planck data. This result is con- 
firmed by other papers of this series. Planck in combination 
with WMAP 9-year large angular scale polarization (WP) yields 
D.K - -0.006 ± 0.018 at 95% CL by combining temperature and 
lensing information (Planck Collaboration XVI, 2013; Planck 
Collaboration XVll, 2013). The bispectral non-Gaussianity pa- 
rameter /nl measured by Planck is consistent with zero (Planck 
Collaboration XXIV, 2013). These results are consistent with 
zero spatial curvature and a small value of /nl, as predicted in 
the simplest slow-roll inflationary models. 

A key Planck result is the measurement of the scalar pertur- 
bation spectral index. Planck+WP data give n^ - 0.9603 ±0.0073 
(and n, = 0.9629 ±0.0057 when also combined with BAG). This 
result disfavours the Harrison-Zeldovich (HZ) n^ = I model at 
more than 5cr. Even in extended cosmological models, the HZ 
spectrum cannot be reconciled with the data. Allowing a general 
reionisation scenario yields ^xlg - 12-5 with respect to ACDM 
for Planck+WP data. When the primordial Helium abundance or 
the effective number of neutrino species are allowed to vary, the 
best-fit of the HZ model to a combination of Planck+WP and 
BAG data is still worse by ^xlg - 4 6, and 8.0, respectively. 

We find no evidence for Planck data preferring a generaliza- 
tion of a simple power law spectrum to include a running of the 
spectral index (dws/dlnA: = -0.0134 ± 0.0090) or a running of 
the running {(fnjdlnk^ = 0.020+°°]^ with Planck+WP). In a 
model admitting tensor fluctuations, the 95% CL bound on the 
tensor-to-scalar ratio is ro.oo2 < 0.12 (< 0.11) using Planck+WP 
(plus high-^). This bound on r implies an upper limit for the in- 
flation energy scale of 1.9 x 10'^ GeV or, equivalently, for the 
Hubble parameter < 3.7 x 10"^ Mpi, at 95% CL. 

The degeneracy between n^ and r, which plagued previous 
CMB measurements, is now removed by the Planck precision 
in the determination of the highest acoustic peaks. Inflaton po- 
tentials with a concave shape are favoured and occupy most of 
the 95% confidence region allowed by Planck+WP in the n^- r 
plane. Models with an exponential potential, a monomial poten- 
tial with a power larger than two, or hybrid models driven by 
a quadratic term are disfavoured at more than 95% confidence. 
The quadratic large field model, in the past often cited as the 
simplest inflationary model, is now at the edge of the 95% CL 
contours aUowed by Planck+WP+\\ig\\-{ CMB data. 

A Bayesian parameter estimation and model comparison 
analysis of a representative sample of single-field slow -roll mod- 
els shows that Planck is able to discriminate between these mod- 
els with results that are robust even when a broad set of entropy 
generation scenarios are allowed. In addition to confirming the 
exclusion of the potential, the Bayesian evidence computed 
from the Planck data provides significant odds (logarithms of the 



Bayes factor of 5 or lower relative to ACDM) against large- 
field models compatible with previous cosmological data, such 
as the (p^ potential, and two-parameter potentials such as natural 
inflation and the hilltop potential. As presented in Sect. 5, Planck 
establishes strong constraints on the parameter values of specific 
inflationary scenarios. For example, the scale parameter of the 
natural inflation potential is constrained to be log(//Mpi) > 1.1 
(95% CL), improving upon the WMAPl limit on / by a factor 
of two. The Planck data limit the possibilities for the unexplored 
physics between the end of inflation and the beginning of the 
radiation-dominated era. Data-driven constraints are obtained on 
Wint, the effective equation of state in the post-inflationary era. 
Particularly for the disfavoured models listed above, their pa- 
rameters are pushed to unnatural values Wint > 1/3 in order to 
become more compatible with the data. 

Using an essentially exact numerical calculation of the pre- 
dicted primordial spectrum, we reconstruct the observable win- 
dow of the inflaton potential, expanding the potential as a Taylor 
series up to a fixed order. For an observable potential described 
by a polynomial of order three, the reconstruction agrees well 
with the slow-roll predictions. If a quartic term is allowed, the 
result deviate from the slow-roll prediction, because the Planck 
data favour a slightly smaller amplitude for the Sachs-Wolfe 
plateau relative to the ^ > 40 part of the power spectrum than 
the best-fitting minimal ACDM model with a power law primor- 
dial spectrum. A potential with an fourth order polynomial can 
fit this feature, thus reducing the effective x\fi by approximately 
four 

A penalized-likelihood reconstruction of the primordial 
power spectrum shows hints of structure at modest estimated sta- 
tistical significance. Possible sources of unaccounted errors are 
discussed. Parameterized models producing superimposed oscil- 
lations (possibly motivated by deviations from the Bunch-Davies 
vacuum state, axion monodromy, or a sharp step in the inflaton 
potential) improve the^^^j. by roughly 10, where three extra pa- 
rameters have been added. However, a Bayesian model compari- 
son analysis does not strongly favour the model with oscillations 
over the standard featureless power spectrum. With Planck po- 
larization data, a more conclusive result on superimposed oscil- 
lations is expected. 

We combine power-spectrum constraints with those on the 
non-linearity parameter /jl (Planck Collaboration XXIV, 2013) 
to constrain single field inflation with generalized Lagrangians, 
in which non-Gaussianities are larger than those predicted by the 
simplest slow-roll inflationary models. We show how the lim- 
its on the inflation sound speed derived in Planck Collaboration 
XXIV (2013) are crucial to constrain slow-roll parameters for 
generalized Lagrangians. We also show how particular exam- 
ples of DBl Inflation and fc-inflation can be constrained by this 
combination of Planck data. 

We test the hypothesis that the primordial cosmological 
perturbations were exclusively adiabatic. We analyse all non- 
singular (i.e., non-decaying) isocurvature modes arbitrarily cor- 
related to the adiabatic mode, using a parameterization where the 
isocurvature contributions are specified at two scales. The oscil- 
latory pattern in the Planck temperature spectrum is compatible 
with purely adiabatic perturbations, and therefore constrains any 
isocurvature contribution to be small at those multipoles. As a 
consequence, axion and curvaton scenarios, in which the CDM 
isocurvature mode is uncorrected or fully correlated with the 
adiabatic mode, are not favoured by Planck. The upper bounds 
on the isocurvature fraction aik - 0.05 Mpc"' are 0.039 for the 
axion, and 0.0025 for the curvaton, at 95% CL. However gen- 
eral models with an arbitrarily correlated mixture of adiabatic 
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and (CDM or neutrino) isocurvature modes have the freedom to 
lower the Sachs-Wolfe plateau relative to the high-/" spectrum, 
and reduce the effective xlff by more than four. 

The simplest inflationary models have passed an exact- 
ing test with the Planck data. The full mission data includ- 
ing Planck's polarization measurements will help answer fur- 
ther fundamental questions, including the possibilities for non- 
smooth power spectra, the energy scale of inflation, and exten- 
sions to more complex models. 
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Appendix A: Sampling tlie Hubbie flow functions 

In this appendix we briefly review how to constrain slow-roll in- 
flation by sampling the Hubble flow functions (HFFs), and dis- 
cuss how well the results agree with those derived by sampling 
directly the parameters InAs, n^, r, and dws/dln^. This method 
fully exploits an analytic perturbative expansion in terms of the 
HFF for the primordial spectra of cosmological fluctuations dur- 
ing slow-roll inflation (Stewart & Lyth, 1993; Gong & Stewart, 
2001; Leach et al., 2002), which self-consistently extends to 
highest order the first terms presented in Eqs. 13-19. Since z" /z 
in Eq. 7 and a" /a in Eq. 9 can be rewritten exactly in terms of 
the Hubble flow functions, the HFF hierarchy, rather than the po- 
tential one, is best suited for this purpose. The slow-roll analytic 
power spectra have been calculated up to second order using the 
Green's function method (Gong & Stewart, 2001; Leach et al., 
2002). Other approximations are available in the literature, in- 
cluding WKB (Martin & Schwarz, 2003), uniform approxima- 
tion (Habib et al., 2002), or the method of comparison equations 
(Casadio et al., 2006). 

The dependence of the amplitudes in Eqs. 13 and 14 in terms 
of HFF is given by 

Ax^Axoe'^o, (A.l) 
where X = s,t, and b^o , bto are: 

Z7sO = - 2 (C -H 1) ei - Ce2 + (-2C + y - 7) 

+ (f - l) 4 + (-^' - 3X + ^ - 7 + A,o) €162 (A.2) 

+ + fi +A,o)e2e3, 

^o=-2(C+l)ei+(-2C+f -7)e2 

+ (-C2-2C+f5-2 + A,o) 6162, 

and C = In 2-i-yE-2 x -0.7296 (yn is the Euler-Mascheroni con- 
stant). At second order the coefficients of the expansion depend 
on the approximation scheme: X - C and Aso = Ato = apply 
for the Green's function method (Gong & Stewart, 2001; Leach 
et al., 2002); X = D = l/3-ln3, A^o = (Z)-C)(D-Hln2)- 1/18, 
Ato - 2D(D - C) - 1/9 apply for for the method of comparison 
equations (Casadio et al., 2006). As predicted by the consistency 
relation, A;,o = IGeiAno- 

The full perturbative expressions up to second order in HFF 
for the spectral indices and the running of the indices are 

- 1 = -2ei - 62 - 2ef - (2 C -H 3) e, 62 - Cez^a, (A.4) 



(A.3) 



dn^/dlnk - -2e\e2 - £2^^,, 

til - -2ei - 2e\ - 2 {C + I) ei 62, 
dnt/dlnA; - -2e\e2. 



(A.5) 
(A.6) 
(A.7) 



We now constrain these parameters using the Planck+WP 
data. If we restrict ourselves to first order, we obtain ei < 0.0074 
at 95% CL and £2 = 0.030!°{;i° at 68% CL. At second or- 
der with the Green's function method, we obtain ei < 0.013 
at 95% CL, 62 = 0.043!|J;JJ|;^ and 63 = 0.36![;;^^ at 68% CL. 
The comparison of results at first and second order is shown in 
Fig. A. 1 . Different approximation schemes lead to small differ- 
ences, as Fig. A.2 shows for the Green's function method versus 
the method of comparison equations. Fig. A.3 shows the agree- 
ment between the physical parameters reconstructed from the 
HFF method and those directly sampled as in Sect. 4.3. We con- 
clude that Planck sensitivity has carried the slow -roll inflationary 
analysis to a regime in which differences among analytic predic- 
tions and/or different priors are sub-leading effects up to second 
order. 
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Fig. A.l. Planck constraints on the HFFs (ei , £2, fs) assuming ei- 
ther 63 = and the first-order slow-roll approximation for the 
computation of the primordial spectra, or 63 and the second- 
order slow-roll approximation (HFF GFM). 



Fig.A.2. Comparison of the Planck constraints on the HFFs 
(61,62,63) using the Green's function method (GFM) and the 
method of comparison equations (MCE). 
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Fig.A.3. Planck constraints on the spectral parameters 
Hs ,dns/dlnfe,r. We compare constraints computed with the 
second-order slow-roll approximation, starting from flat priors 
on the HFF parameters at the pivot scale (ei, 62, £3), with those 
obtained directly from , dws/d In A; , r. In the latter case we en- 
force the second-order consistency conditions for the tensor-to- 
scalar ratio and for the running of the tensor spectral index. 



